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Abstract 


CITPUFF,  is  a  puff-type  dispersion  model  that  uses  a  wind  field 
calculated  from  a  complex-terrain  wind  model.  It  accommodates  a 
variety  of  source  types  including  point,  area,  and  line  sources; 
calculates  plume  rise  where  applicable;  and  outputs  a  graphic 
display  of  puff  trajectories  and  concentrations.  The  model  is  com- 
pared against  models  currently  used  for  assessing  air  quality  impacts 
in  complex  topography. 
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MANAGEMENT  IMPLICATIONS 

This  paper  presents  technical  details  of  an  air  pollu- 
tion dispersion  model  implemented  as  a  component  of 
the  Topographic  Air  Pollution  Analysis  System  (TAPAS). 
TAPAS  is  a  user-friendly,  semi-interactive,  collection  of 
computer  programs  implemented  on  the  Colorado  State 
University  CDC-Cyber  computer.  The  programs  are  writ- 
ten in  CDC's  version  of  FORTRAN  IV. 

TAPAS  includes  a  number  of  programs  that  fit  into 
three  categories:  file  management  routines,  models,  or 
input/output  routines.  File  management  routines  allow 
the  user  to  construct  data  files,  model  submit  files,  and 
output  submit  files.  The  input/output  routines  allow  the 
user  to  access  data  and  to  output  graphic  products  to 
various  graphics-capable  hardware.  The  models,  how- 
ever, represent  the  central  value  of  TAPAS.  They  are 
designed  to  simulate  air  flow  in  mountainous  terrain  and 
the  dispersion  of  pollution  released  from  various 
sources  in  mountainous  terrain.  A  primary  objective  of 
TAPAS  is  to  increase  the  utility  of  such  models  by  in- 
cluding them  in  a  user-oriented  environment.  A  second 
objective  is  to  provide  as  wide  an  array  of  models  as 
possible.  This  is  important  because  mountain  air  quality 
and  related  air  resource  problems  are  very  complex.  It 
is  unlikely  that  any  one  approach,  technique,  or  model 
will  be  adequate  for  all  purposes.  Furthermore,  suffi- 
cient data  resources  do  not  exist  to  establish  overall 
validity  of  alternative  techniques.  TAPAS  provides  a 
range  of  techniques  in  a  friendly  environment  so  that  the 
user  can  select  the  most  appropriate  technique  for  a 
particular  application. 

The  CITPUFF  model  described  in  this  paper  is  a 
dispersion  model;  together  with  the  WINDS  model,  it 
provides  an  approximation  of  the  physical  situation 
which  may  be  helpful  in  conducting  certain  types  of 
planning  analyses.  It  is  a  reasonably  economical  system. 
Should  the  user  interest  be  in  more  specific  real-time  or 
operations-oriented  modeling,  the  approach  here  may  be 
useful  but  great  caution  is  appropriate.  The  last  section 
of  this  report  describes  the  limitations  in  the  model 
physics,  which  are  significant,  and  highlights  what  may 
be  appropriate  applications. 

INTRODUCTION 

Predicting  the  concentration  of  air  pollutants  is  a 
problem  that  increases  in  significance  with  industrial 
development  and  population  increases.  In  mountainous 
areas,  particularly  in  the  western  United  States, 
development  is  generally  associated  with  mineral  pro- 
duction and  recreational  opportunities.  The  Clean  Air 


Act  requires  significant  involvement  for  federal  land 
managers  in  the  review  of  prevention  of  significant 
deterioration  (PSD)  permits  (Part  C,  Public  Law  95-95). 
In  addition,  federal  land  managers  are  responsible  for 
managing  pollution  their  activities  generate,  like  smoke 
from  prescribed  burning.  For  both  of  these  air  quality 
models  capable  of  predicting  the  concentration  of  pollu- 
tion that  would  result  from  a  source  in  mountainous  ter- 
rain are  needed.  A  gaussian  puff  model,  the  Complex 
Interactions  with  Terrain-Puff  Model  (CITPUFF),  was 
designed  to  partially  meet  this  need. 

To  predict  concentrations  in  complex  topography,  it  is 
necessary  to  estimate  (1)  the  meteorology  of  the  region  in 
some  detail;  (2)  the  diffusion  of  pollution;  and  (3)  the 
removal  and/or  chemical  transformation  of  pollutants. 
To  support  these  estimates,  the  Topographic  Air  Pollu- 
tion Analysis  System  (TAPAS)  has  been  developed  (Fox 
and  Fosberg  1977,  Childs  and  Marlatt  1981,  Fox  et  al. 
1984).  TAPAS  includes  a  series  of  air  quality  related 
models  and  procedures.  Because  applications  are  an- 
ticipated in  data-poor  areas,  wind  models  are  among  the 
most  significant  elements  of  TAPAS.  WINDS  (Fosberg  et 
al.  1976)  is  a  TAPAS  model  that  is  capable  of  providing  a 
grid-based,  diagnostic,  two-dimensional  distribution  of 
wind  speed  and  direction.  WINDS  provides  an  analysis 
of  the  distribution  of  winds  in  complex  terrain  based  on 
an  assumed  (or  measured)  primary  influencing  wind. 
The  WINDS  model  generated  wind  field  can  be  used  as 
the  basis  for  estimating  plume  trajectories  in  the  pollu- 
tion transport  and  dispersion  modules  of  TAPAS.  This 
document  describes  the  primary  TAPAS  pollution  trans- 
port and  dispersion  module,  CITPUFF. 

CITPUFF  is  a  Lagrangian  "puff-type"  diffusion 
model  based  on  the  Start  and  Wendell  (1974)  model, 
MESODIF  (also  Bass  et  al.  1979).  Pollution  released 
from  a  source  is  transported  by  the  wind  while  it  dif- 
fuses as  a  result  of  the  atmospheric  turbulence.  The 
details  of  the  source  must  be  described  in  terms  of  the 
source  location,  type  of  pollution,  rate  of  emission,  and 
the  physical  source  parameters  (point  source,  area 
source,  height  of  stacks,  etc.).  CITPUFF  calculates  the 
appropriate  plume  rise  using  EPA-recommended  proce- 
dures (EPA  1978,  1980).  An  initial  source  dispersion 
parameter  is  used  as  an  option  to  accommodate  a  num- 
ber of  different  source  types,  ranging  from  an  isolated 
point  source,  to  a  time-dependent  forest  burn,  to  an 
area  source  such  as  a  surface  mine.  The  puff  is  dif- 
fused according  to  conventional  theories  of  atmo- 
spheric turbulence.  Empirical  classes  of  atmospheric 
stability  used  in  the  model  are  related  to  the  level  of  at- 
mospheric turbulence.  The  user  may  select  from  three 
options  for  dispersion,  including  a  fixed  22.5°  horizon- 
tal spread,  standard  Pasquill-Gifford-Turner  disper- 
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sion  coefficients,  and  a  convective  boundary  layer 
based  scheme  using  altered  Pasquill-Gifford  coeffi- 
cients. The  CITPUFF  model  also  allows  a  linear  deple- 
tion of  mass  from  the  puff.  Puffs  are  released  at  a  user- 
selected  rate  and  tracked  until  they  are  transported  off 
the  wind  field  grid.  Concentration  is  accumulated  at 
receptors  located  at  each  computational  grid  point  to 
determine  the  ground-level  concentration  of  pollutants. 

First,  the  document  describes  the  general  theory  on 
which  CITPUFF  is  based.  Then  it  develops  the  details  of 
CITPUFF,  including  the  grid  system,  meteorology,  emis- 
sions configurations,  terrain  height  adjustment,  diffu- 
sion coefficients,  mixing  depth  adjustments,  and  puff 
sampling  functions.  The  last  section  identifies  how 
CITPUFF  performs  relative  to  other  models  and  model- 
ing techniques  and  discusses  appropriate  limitations 
on  model  use. 


CITPUFF  MODEL— GENERAL  THEORY 

When  the  principle  of  conservation  of  mass  is  ap- 
plied in  the  atmosphere,  together  with  a  Fickian  ap- 
proximation for  molecular  diffusion,  it  yields  the 
following  form  of  the  diffusion  equation: 

^  +V(XV)  =  V(DVX)  +  ES  [1] 

where 

x(x,t)  =  instantaneous,  realized  concentration  at 
time  t  at  position  x  relative  to  a  given 
origin, 

x-(x,y,z)  =  Cartesian  coordinates, 
V  =  velocity  field, 
D  =  molecular  diffusivity, 
ES  =  sum  of  sources  and  sinks. 

In  the  present  case,  ES  is  replaced  by  the  term  -kx 
which  may  be  used  to  account  for  such  nonconservative 
effects  as  first-order  removal,  conversion  of  gaseous  to 
particulate  pollutants,  and  radioactive  decay.  In  such 
cases  the  rate  of  loss  of  the  diffusing  material  is  as- 
sumed proportional  to  its  concentration,  the  factor  of 
proportionality,  k,  being  a  "decay  constant"  of  dimen- 
sions s_1. 

In  turbulent  flow  both  fluid  velocity  and  tracer  con- 
centration are  random  variables.  These  variables  may 
be  separated  into  ensemble  mean  and  fluctuating  com- 
ponents. To  focus  attention  on  the  effects  of  complex 
topography  on  dispersion,  the  approach  of  Fosberg  and 
Fox  (1978)  is  extended  and  the  fluctuating  components 
are  considered  to  be  composed  of  topographically  in- 
duced effects  and  the  usual  flat-terrain  turbulent  fluc- 
tuation, as  follows: 

V  =  V+  V  +  V 

x  =  x  +  x  +  x' 


where 

V,x  =  ensemble  mean  values, 

V,x  =  terrain-induced  fluctuations  about  the 
mean, 

V,'x'  =the  remaining  small-scale  variations  from 
the  ensemble  mean,  generally  considered 
as  flat-terrain  turbulent  fluctuations. 

Taking  ensemble  averages  of  both  sides  of  equation  [1] 
yields 

^  +V(xV)  +  V(x~V)  +  V(x'V')  =  V(DVx)-kx  [2] 


The  velocity-concentration  correlation  x'V'  represents 
flux  by  turbulent  movements  and  can  be  related  to  mean 
flow  variables  with  a  Fickian  diffusion  approximation. 
Eddy  diffusivities,  K  =  (K^,  Kz),  can  be  substituted  into 
equation  [2] 

^  +Wx  =  V(KVx)-kx-V(xV)  [3] 

where  molecular  diffusion  has  been  neglected  relative  to 
turbulent  diffusion,  and  the  second  term  on  the  left-hand 
side  of  equation  [2]  has  been  simplified  by  using  the 
ensemble  mean  form  of  the  continuity  equation  for  an 
incompressible  fluid  (VV  =  0). 

The  concentration-velocity  correlation,  xV,  represents 
flux  by  terrain-induced  motion.  In  order  to  "directly" 
close  equation  [3]  one  can  postulate  a  Fickian-type  diffu- 
sion approximation  to  relate  x  V  mean  flow  variables, 
with  Cartesian  components  of  the  form: 

XU=  -C  — ±-  ,  Xv=  -Cy~rf-  ,  XW=  -Cz  — — 
ox  y   oy  oz 

The  terrain-induced  "apparent  diffusivities"  (Cx,  Cy,  Cz) 
will  depend  both  on  nonturbulent,  organized,  circulation 
patterns  induced  by  the  terrain,  and  on  topographically 
induced  turbulent  enhancement.  In  both  cases,  there  will 
be  a  direct,  or  at  least  implicit,  dependence  on  differential 
thermal  convection  from  mountain  slope  surfaces,  coupl- 
ed with  mechanical  interactions  of  the  terrain  with  the 
mean  winds. 

The  nonturbulent  component  accounts  for  terrain- 
induced  kinematic  constraints  of  the  flow,  and  their 
subsequent  effects  on  horizontal  and  vertical  dispersion. 
These  kinematic  constraints  are  produced  by  distortions 
caused  by  air  flow  around  terrain  obstacles  that  produce 
phases  of  acceleration  and  deceleration  that  can  be 
regarded  as  independent  of  turbulent  diffusion  to  the  first 
order. 

Terrain-induced  kinematic  constraints  on  plume 
dispersion  and  plume  trajectories  have  been  examined  by 
Hunt  and  Mulhearn  (1973),  Egan  (1975),  and  Hunt  et  al. 
(1978)  using  the  theory  of  turbulent  plumes  embedded 
within  potential  flow  fields.  Bass  et  al.  (1981)  extended 
this  approach  to  slightly  stable  flows  by  an  empirical  ap- 
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proximation  derived  from  stratified  tow-tank  experiments 
(Snyder  1981),  as  well  as  for  objects  of  arbitrary  crosswind 
aspect  ratio.  Bass  et  al.  (1981)  present  a  comprehensive 
review  of  the  work  performed  by  using  this  approach,  as 
well  as  a  discussion  of  the  associated  limitations.  In  par- 
ticular, this  report  provides  a  detailed  description  of  the 
circumstances  and  limitations  associated  with  a 
simplistic  inclusion  of  flow  fluid  distortions  on  a  Gaus- 
sian plume  model. 

Recently,  a  considerable  amount  of  research  has  been 
directed  toward  understanding  the  behavior  of  plumes  in 
simplified  topographic  settings.  Laboratory  simulation 
and  actual  field  studies  have  been  compiled  for  Cinder 
Cone  Butte  in  Idaho,  a  small,  Gaussian-shaped,  100-m- 
high  hill.  The  reports  from  this  study  (Lavery  et  al.  1982, 
Lavery  et  al.  1983,  Strimaitis  et  al.  1983)  suggest  that 
dispersion  coefficients  for  use  in  modified  Gaussian 
models  should  be  calculated  from  measurements  of  the 
local  turbulence  level.  This  is  not  a  particularly  surprising 
result,  and  it  suggests  that  emphasis  should  be  placed  on 
proper  measurement,  or  simulation  of  the  wind  field  and 
its  associated  turbulence.  These  studies  have  also  clearly 
established  the  concept  of  the  dividing  streamline  in 
stable  flows.  This  is  a  level  in  the  atmosphere  where  the 
kinetic  energy  an  air  parcel  contains  is  matched  by  the 
potential  energy  required  to  lift  the  parcel  against  the  den- 
sity gradient.  This  height,  in  a  simple  situation  with 
uniform  wind  and  temperature  profiles,  is  given  by 

hd  =  (l-Fr)h 


F  =  -iL- 
r  Nh 


N 


-\T  ~dz~J 


where 

hd  =  dividing  streamline  height, 
u  =  uniform  velocity  upstream  of  the  hill, 
h  =  height  of  the  topography, 
d  =  potential  temperature, 
g  =  acceleration  of  gravity. 

This  height  separates  two  vastly  different  flow 
regimes.  Above  this  level,  air  will  flow  over  the 
topography;  below  this  level,  air  essentially  stagnates. 
This  research  clearly  demonstrates  the  need  to  make 
local  measurements  in  complex  terrain.  It  also  shows 
that  even  simple  situations  in  complex  topography  lead 
to  remarkable  complexity  in  the  flow.  As  a  final  objec- 
tive, this  EPA-supported  research  is  expected  to  yield 
improved  formulations  for  dispersion  coefficients  and 
stability  effects  in  complex  topography.  In  view  of  this, 
CITPUFF  employs  a  modular  approach  allowing  the  in- 
corporation of  new  procedures  as  they  become 
available. 

Within  the  context  of  a  puff  model,  terrain-induced 
alterations  of  the  puff/plume  centerline  can  be  handled 
via  the  wind  field  input  that  is  used  to  transport  the  puff. 
The  effect  of  terrain  on  dispersion  by  kinematic  effects 
and  the  effects  of  local  turbulence  can  be  incorporated 


by  adjusting  the  flat-terrain  dispersion  coefficients.  The 
latter  adjustments  could  perhaps  be  achieved  by 
postulating  relationships  between  the  diffusivities  Cx,  C  , 
and  C2  and  variables  such  as  streamline  convergence 
and  divergence,  which  characterize  phases  of  accelera- 
tion and  deceleration.  This  approach  must  be  merged 
with  an  understanding  of  turbulence  generation  via 
these  mechanisms  and  is  beyond  the  scope  of  this  report. 
Unfortunately,  the  current  state  of  knowledge  in  this 
area  (Egan  1984)4  is  inadequate.  The  current  discussion 
is  included  simply  to  illustrate  that  a  theoretical  basis 
and  a  mechanism  for  incorporating  these  effects  into  a 
Gaussian-type  model  by  the  dispersion  relations  exist. 
The  discussion  that  follows,  however,  will  use  the  flat- 
terrain  dispersion  relations  and  account  for  terrain  ef- 
fects via  the  wind  field  and  a  terrain  height  adjustment 
algorithm  because  the  state  of  knowledge  at  present 
does  not  support  alteration  of  dispersion  coefficients 
(Strimaitis  et  al.  1983). 
Equation  [3]  can  now  be  written  in  the  form: 


dx 

at 


+  u- 


3x 


d 

dx 


dx 
dx 


d 

ay 


K 


,  dX 
y  dy 


where 


V  =  (u,0,0) 

The  final  term  on  both  sides  of  equation  [4]  can  be 
eliminated  by  the  change  of  variable 

c  =  x  exp  ( -  kx/u)  exp  (xu/2Kx0  exp  ( -  \itl4K$ 
to  obtain 


JLiL 

at 


dx 


K  '4^1 
y  dy) 


[5] 


A  solution  to  equation  [5]  for  an  instantaneous  point 
source  of  strength  Q  (units  of  mass),  generated  at  time 
t  =  0  located  at  the  origin  of  the  Cartesian  coordinate 
system  was  first  obtained  by  Roberts  (1923)  as 


(27r)3'2axaa 


exp 


xz  yz  zz 
2a2      2a2  2a2 


[6] 


where  it  has  been  assumed  that  a  Gaussian  distribution 
in  x,  y,  and  z  for  the  concentration  exists,  and  such  a 
distribution  is  related  to  the  "eddy"  diffusivity  (assumed 
constant)  by 


ax  =  2tKx  =  2t(Kx  +  Cx) 
a2  =  2tKy'  =  2t(Ky  +  Cy) 
<T*=2tK'=2t(K+CJ 


[7] 


'Egan,  B.  A.  1984.  Workshop  report  on  air  quality  modeling  in 
complex  terrain— Keystone,  Colorado,  May  1983.  American  Mete- 
orological Society,  Boston,  Mass.  (In  press). 
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The  solution  is 


X  = 


(2irf'\ayaz 


exp 


_Ut]2    |     y2     |     Z2  ' 
2ff? 


2a? 


2a2 


exp(-kx/u) 


[8] 


for  the  case  of  a  puff  from  an  instantaneous  point  source 
that  is_advected  in  the  x  direction  with  a  mean  wind 
speed  u,  while  losing  mass  at  an  exponential  rate  k  as  a 
result  of  the  assumed  first-order  decay. 

For  a  puff  generated  from  a  point  source  with  an  effec- 
tive emission  height  h  above  ground  level,  equation  [8] 
becomes 


X(x,y,z)=       ,  exp 

>         V  Z 

+  exp 


y  z 

-(z  +  h)2 


2a2 


exp 


(z-h)2 


2a? 


2a2 


[9] 


exp  ( -  kx/u) 


where  r2  =  (x  —  ut)2  +  y2,  and  where  horizontal  Gaussian 
symmetry  has  been  assumed  (i.e.,  ax  =  ay).  In  this  form  the 
actual  source  (with  stack  height  hs)  is  located  at  the  point 
(0,0,  hs).  The  second  term  in  the  large  bracket  results  from 
an  image  source  (of  the  effective  source)  at  (0,0,  —  hj  to 
represent  the  effect  of  the  ground. 

The  resulting  ground-level  concentration  for  a  discrete 
puff,  in  the  form 


2Q 


x(x,y.0)=  (27r)3;2ff2g  exp 


2a2 


exp 


2a2 


exp 


kx 


[10] 


is  used  as  the  basic  building  block  for  the  puff  compo- 
nent of  the  model. 

Source  emissions  are  simulated  as  a  time  series  of 
discrete  puffs.  At  each  time  step,  the  puff  centroids  are 
advected  by  a  two-dimensional  Lagrangian  trajectory 
function.  The  standard  deviations  of  the  assumed  Gaus- 
sian concentration  distributions  are  taken  to  be  power 
law  functions  of  travel  distance/time. 

The  extension  to  multiple  sources  is  based  on  the 
assumption  that  the  total  concentration  field  of  a  pollu- 
tant can  be  obtained  by  superposition  of  the  concentra- 
tion fields  for  each  individual  source. 

Line,  area,  and  volume  source  distributions  can  be 
simulated  by  the  use  of  virtual  point  sources  and  the 
specifications  of  initial  values  of  the  Gaussian  standard 
deviations  (ayo  and  azo)  characteristic  of  the  particular 
source  distribution.  Care  must  be  exercised  in  modeling 
cases  where  these  dimensions  may  vary  with  temporal 
or  local  spatial  changes  in  wind  direction  and  speed. 
When  local  spatial  changes  are  evident,  it  may  be 
necessary  to  represent  a  large  area  source  by  more  than 
one  puff  in  order  to  properly  represent  the  flow,  since 
the  trajectory  function  only  considers  the  centroid  of 
each  puff. 


The  CITPUFF  model  is  designed  to  predict  pollutant 
concentrations  at  specific  receptors.  At  any  such  recep- 
tor the  model  calculates  the  integrated  concentration, 


[u] 


where  T  refers  to  the  dosage  time  and  ^s  is  often  refer- 
red to  as  the  dosage  or  exposure,  and  the  subscript  s  is 
introduced  as  a  reminder  that  all  sources  and  puffs  are 
included  in  the  calculation.  The  average  concentration 
over  the  time  interval  T  is  then  found  by  dividing  ^s  by  T. 
The  dosage  integral  in  equation  [11]  is  approximated  in 
CITPUFF  by  the  use  of  a  sampling  function.  The  method 
of  calculation  described  above  will  also  simulate  the 
concentration  field  from  a  continuous  source,  provided 
the  series  of  discrete  puffs  overlap  along  the  variable 
plume  axis  with  sufficient  density.  The  last  section  con- 
siders simulation  of  the  continuous  plume  in  more  detail, 
and  also  compares  the  predictions  of  CITPUFF  to  those 
of  other  models. 


CITPUFF  MODEL— SPECIFIC  MODEL 
COMPONENTS 

Grid  System  and  Basic  Computational  Unit 

The  CITPUFF  model  uses  a  Cartesian  coordinate  sys- 
tem for  both  input  and  output.  The  maximum  allowable 
size  of  the  computational  grid  in  the  current  version  of 
CITPUFF  (as  implemented  on  the  CSU  Cyber  720  com- 
puter) is  approximately  6400  grid  points  (e.g.,  an  80  by 
80  grid  area). 

The  wind  field  required  to  drive  the  CITPUFF  model 
must  be  specified  in  the  form  of  horizontal  (u  and  v)  wind 
vector  components  at  each  Cartesian  grid  point  of  the 
computational  grid.  The  wind  field  can  be  updated  after 
each  basic  computational  unit  (1-hr  intervals  in  the  cur- 
rent version  of  CITPUFF).  The  source  of  the  wind  fields 
used  in  CITPUFF  and  the  subsequent  interpolation 
scheme  used  for  advecting  puffs  are  outlined  in  the  next 
section. 

The  basic  computational  unit  in  the  current  version  of 
CITPUFF  is  1  hr.  The  model  is  organized  so  that  the 
minimum  time  average  that  can  be  obtained  directly  is  a 
1-hr  average.  Emission  characteristics  and  meteorology 
can  be  updated  every  hour,  with  the  input  values  taken 
as  representative  of  1-hr  averages  for  the  given  period 
of  time. 

Factors  such  as  puff  release  rate,  sampling  frequen- 
cy, time  and  length  of  basic  advection  steps,  are  de- 
signed to  support  the  1-hr  basic  computational  unit.  The 
CITPUFF  model  is  flexible,  and  if  required  can  be  ad- 
justed for  alternative  basic  computational  units. 


Meteorology 

The  CITPUFF  model  requires  the  following  meteoro- 
logical inputs: 
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a.  Stability  class  and  mixing  depth  specified  at  hourly 
intervals 

b.  Horizontal  (u,v)  wind  components  (at  a  standard 
height  above  ground  level)  specified  at  each  of  the 
wind  field  grid  points  at  hourly  intervals 

c.  Ambient  temperature  (in  °K)  to  calculate  plume 
rise  for  any  sources  in  the  analysis  area  at  hourly 
intervals. 

In  the  current  version  of  CITPUFF  a  single  stability 
class  mixing  depth  and  ambient  temperature  are  spe- 
cified for  the  analysis  area  being  considered.  The 
gridded  wind  components  (u  and  v)  currently  used  to 
drive  CITPUFF  are  produced  by  the  complex-terrain 
wind  field  model  WINDS  (Fosberg  et  al.  1976).  It  is  also 
possible  to  use  the  u  and  v  components  from  a  three- 
dimensional  wind  model,  which  like  WINDS  is  a  compo- 
nent of  the  Topographic  Air  Resources  Analysis  System 
(TAPAS)  (Fox  et  al.  1984).  The  WINDS  model  is  a  grid- 
based,  two-dimensional  diagnostic  model.  The  model 
relies  on  simplified  solutions  to  complex  equations  that 
characterize  the  major  controlling  physical  processes 
that  influence  wind  flow  over  complex  terrain.  The 
model  does  not  attempt  to  solve  time-dependent,  predic- 
tive, atmospheric  fluid  and  thermodynamic  equations 
that  require  large  initial  data  fields  and  extensive  com- 
puter resources.  Instead,  the  model  simulates  the  two- 
dimensional,  horizontal  wind  flow  over  complex  terrain 
by  combining  the  terrain-induced  flow  concept  de- 
scribed by  Anderson  (1971)  and  the  thermally  induced 
wind  flow  processes  described  by  Fosberg  et  al.  (1972). 
The  WINDS  model  handles  subsynoptic  situations 
where  a  characteristic  prevailing  wind  can  be  used  to 
initialize  the  entire  grid.  Intermittent,  decoupled,  or 
separated  flow  simulations  are  not  within  the  scope  of 
the  model.  Hence,  there  are  limitations  on  the  severity  of 
topographic  variation  that  can  be  addressed  by  the 
model.  The  model  was  designed,  and  is  best  applied,  to 
accommodate  grid  point  spacing  of  one-half  to  several 
kilometers. 

The  primary  model  inputs  are: 

•  Elevation  at  each  grid  point 

•  Vegetation/land  cover  class  at  each  grid  point 

•  Temperatures  and  associated  heights  of  two  layers 
of  the  atmosphere  that  define  the  atmospheric  tem- 
perature profile 

•  Influencing  wind  direction  and  wind  speed  (and 
any  available  wind  observations). 

The  WINDS  model  calculates  a  resultant  wind  field  in 
two  phases.  Phase  1  accounts  for  the  effect  of  major  ter- 
rain features  on  the  influencing  wind  throughout  the 
depth  of  the  boundary  layer.  Phase  2  modifies  the  wind 
field  calculated  in  Phase  1  by  accounting  for  wind  per- 
turbations caused  by  the  actual  terrain,  thermal,  and 
frictional  forces  within  a  shallow,  near-surface  air 
layer.  In  both  phases,  mass  is  conserved  within  defined 
horizontal  and  vertical  dimensions  and  nonlinear  forces 
are  not  calculated. 

The  primary  output  of  the  model  used  in  CITPUFF  is  a 
digital  file  that  contains  both  the  u  and  v  wind  com- 
ponents at  each  grid  point  of  the  computation  grid  for  a 


selected  analysis  area.  The  WINDS  model  as  well  as 
other  TAPAS  models  outputs  are  formatted  in  TAPAS- 
structured  files  (Childs  and  Marlatt  1981).  CITPUFF  is 
designed  to  directly  read  and  interpret  TAPAS- 
structured  files. 

Trajectory  Calculation 

The  trajectory  calculation  requires  two  gridded  wind 
field  data  sets  that  represent  conditions  occurring  1  hr 
apart.  These  wind  component  data  sets  are  used  to  con- 
struct trajectories  at  a  particular  place  on  the  grid  at  a 
particular  time,  by  performing  a  bilinear  interpolation 
in  space,  followed  by  a  linear  interpolation  in  time. 

The  trajectory  for  the  u  component  is  illustrated 
following  the  original  presentation  in  Hales  et  al.  (1977). 
Let  the  times  of  two  consecutive  wind  fields  be  t  and 
t  +  Atw,  respectively.  The  interval  between  wind  maps, 
Atw,  is  1  hr  in  the  current  version  of  CITPUFF.  Referring 
to  figure  1,  it  is  necessary  to  calculate  u(t  +  At,p,q),  given 
u(t,0,0),  u(t,l,l),  u(t  +  Atw,l,l),  u(t  +  Atw,0,0),  etc.,  where  At 
is  the  advection  time  interval  (time  step)  and  p  and  q  rep- 
resent grid  points.  In  this  case, 


u(t  +  At.p.q): 


At 


At 


[(l-p)(l-q)u(t,0,0) 


+  (1  -  p)qu(t,0,l)  +  p(l  -  q)u(t,l,0)  +  pqu(t,l,l)] 

+  "  P)(1 "  ^  +  Atw'°'0)  +  (1 "  P)Qu(t  +  Atw0,1) 


+  p(l  -  q)u(t  +  Atw,l,0)  +  pqu(t  +  Atw,l,l)] 

The  same  procedure  is  used  for  determining  v(t  +  At,p,qj. 

Winds  computed  in  this  manner  from  grid  values  are 
used  in  a  two-step  iteration  method  to  calculate  trajec- 
tory steps.  The  manner  in  which  the  trajectory  step  is 
calculated  from  time  t  to  time  t  +  At  for  a  particle  located 
at  [x(t),y(t)]  in  figure  1  is  such  that  a  particle  so 
transported  would  arrive  at  (x^yj  at  time  t  +  At  if  the 
wind  remained  constant.  However,  because  the  wind 
would  change  along  the  trajectory,  the  actual  location 
[x(t  -I-  At),y(t  +  At)]  of  the  particle  would  be  at  some  nearby 
point.  To  better  approximate  this  location,  the  wind  at 
(Xi.yJ  at  time  t  +  At  is  also  calculated  and  used  to 
calculate  a  second  advection  step  to  (x2,y2).  The 
[x(t  +  At),y(t  +  At)]  is  assumed  to  lie  halfway  between 
[x(t),y(t)]  and  (x2,y2).  The  equations  for  calculating 
x(t  +  At)  are  as  follows: 

Xl  =  x(t)  +  u(t,x(t),y(t))At 

x2  =  xt  +  u(t  +  At.Xj.yJAt 

x(t  +  At)  =  0.5[x(t)  +  x2] 

In  exactly  the  same  way,  y(t  +  At)  is  calculated  using  the 
v  component  of  the  wind. 


Puff  Dispersion 

The  CITPUFF  model  contains  three  user  options  for 
calculating  puff  dispersion:  Pasquill  stability  (Hales  et 
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al.  1977),  Power  Plant  Siting  Program  (PPSP)  dispersion 
coefficients  (Brower  1982),  and  22.5°  sector  dispersion 
(EPA  1981). 


Pasquill-Gifford  Stability  Based  Dispersion 

The  Pasquill-Gifford-Turner  stability  based  dispersion 
approach  is  the  primary  method  used  in  most  CITPUFF 
applications  (Hanna  et  al.  1982).  The  diffusion  of  the 
puffs  relative  to  their  trajectory  is  described  by 
distance-dependent  values  of  oy  and  ax  for  distances  of  up 
to  approximately  100  km.  The  growth  of  a  puff  along  the 
puff  trajectory  distance  s  from  the  source  is  represented 
by: 

dayi 

ay(s  +  As)  =  ay(s)  +  As— |s=As/2 

az(s  +  As)  =  az(s)  +  As—  |s=As/2 

Illustration  of  wind  component 
interpolation  problem 


(0,1) 

(1,1) 

A 

Y 

u(t  +  At,  p,  q) 

q 

P 

(0,0) 

(1,0) 

Two-step  advection  iteration 

(X,Y,)               (X2,  Y2) 

[x(t  +  At),  Y(t  +  At)] 

[X(t),  Y(t)] 

Figure  1.— Interpolation  of  wind  components  and  analysis  of  the 
two-step  advection  iteration. 


Similar  equations  for  ay  and  oz  as  a  function  of  travel 
time  are  also  available  to  calculate  dispersion  as  a  puff 
travels  beyond  approximately  100  km  from  the  source. 

The  distance  functions  used  in  CITPUFF  for  Pasquill- 
Gifford-Turner  atmospheric  stability  classes  A  (very 
unstable)  through  F  (very  stable)  have  the  following 
form: 

ffy  =  asb   ,  az  =  csd 

where  the  coefficients  a,  b,  c,  and  d  correspond  to  a 
given  Pasquill  stability  category,  as  shown  in  table  1. 
These  coefficients  were  obtained  by  approximating  the 
ay  and  az  dispersion  curves  of  Turner  (1970).  Values  of  s 
and  the  corresponding  values  of  ay  and  az  are  in  meters. 
In  CITPUFF,  these  dispersion  functions  are  implemented 
as  derivatives  along  the  puff  trajectory  path. 

For  travel  distances  greater  than  100  km,  time-based 
formulations  are  more  appropriate  than  distance-based 
formulations  (Hanna  et  al.  1982).  The  CITPUFF  model 
contains  a  time  formulation  given  by  Heffter  (1965): 

<7y  =  0.5t   ,   az  =  V2KJ 

where  t  is  in  seconds,  ay  and  az  are  in  meters,  and  Kz  is  the 
vertical  eddy  diffusivity.  Following  Heffter,  Kz  has  been 
set  at  5  m2/s  in  the  CITPUFF  model.  Selection  of  the  time- 
related  dispersion  parameter  is  a  user  option  in 
CITPUFF.  The  user  assigns  a  time,  depending  on  the 
wind  speeds  in  the  analysis  area,  that  corresponds  to 
the  transport  of  pollutants  approximately  100  km  from 
the  source. 


Power  Plant  Siting  Program  Dispersion 

An  option  exists  in  CITPUFF  to  select  the  dispersion 
coefficients  <ry  and  az  used  in  the  Maryland  Power  Plant 
Siting  Program  (PPSP)  (Brower  1982),  as  adapted  from 
Gifford  (1975).  These  coefficients,  initially  suggested  by 
Briggs,  following,  in  part,  Smith  (1968),  are  based  on  em- 
pirical formulations  derived  from  available  diffusion 
data  from  buoyant  sources  with  ground  and  elevated 
releases.  Briggs  does  not  recommend  their  use  beyond  a 
few  kilometers  from  the  source.  Table  2  presents  the 
downwind  distance-dependent  dispersion  coefficients 
as  a  function  of  stability  as  defined  by  Briggs.  Briggs'  A 
coefficients  are  appropriate  for  most  unstable  at- 
mospheric conditions,  while  Briggs'  F  coefficients  are 
appropriate  for  most  stable  conditions.  To  select  the  ap- 
propriate PPSP  dispersion  parameter,  it  is  necessary  to 
evaluate  which  Briggs  stability  category  is  most  ap- 
propriate for  desired  modeling  conditions.  This  must  be 
done  external  to  the  CITPUFF  model  using  the  method 
described  in  Brower  (1982).  Appendix  A  includes  a  com- 
parison between  PPSP  and  the  Pasquill-Gifford-Turner 
dispersion  approaches. 

The  PPSP  dispersion  parameters,  as  applied  in 
CITPUTT,  have  the  forms  given  in  table  2.  Again  they 
are  implemented  in  differential  forms.  The  user  has  the 
option  of  specifying  a  time  at  which  the  time-dependent 
dispersion  coefficients  are  applied. 
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Table  1  — 

Coefficients  for  Pasquil 

stability  based  dispersion  parameter  formulas 

(from  Hales  et  al.  1977). 

Pasquill 

a 
z 

stability  category 

a 

b 

C 

a 

A  (1) 

0.36 

0.9 

0.00023 

2.10 

B  (2) 

0.25 

0.9 

0.058 

1.09 

C  (3) 

0.19 

0.9 

0.11 

0.91 

D  (4) 

0.13 

0.9 

0.57 

0.58 

E  (5) 

0.096 

0.9 

0.85 

0.47 

F  (6) 

0.063 

0.9 

0.77 

0.42 

Table  2.— Horizontal  (a  and  vertical  (az)  dispersion  coefficients  recommended  by  Briggs  (Gifford 
1975)  as  a  function  of  downwind  distance,  x  (in  meters). 


Sigma  curve 


identifer 

ay 

°z 

1.  BRIGGS  A 

0.22x(1  +0.0001  X)" 

1/2 

0.20x 

2.  BRIGGS  B 

0.16x(1  +0.0001X)" 

1/2 

0.12x 

3.  BRIGGS  C 

0.11x(1  +  0.0001  x)" 

1/2 

0.08x(1+0.0002x)-1'2 

4.  BRIGGS  D 

0.08x(1  +0.0001  X)" 

1/2 

0.06x(1  +0.001 5x)" 1,2 

5.  BRIGGS  E 

0.06x(1  +0.0001X)" 

1/2 

0.03x(1  +0.0003X)-1 

6.  BRIGGS  F 

0.04x(1  +  0.0001  X)" 

1/2 

0.016x(1  +0.003x)-1 

sion.  The  CITPUFF  model  has  incorporated  plume  rise 
formulations  from  the  EPA  Multiple  Point  Gaussian 
Dispersion  Algorithm  with  Optional  Terrain  Adjustment 
(MPTER)  (EPA  1980)  to  calculate  the  effective  emission 
height  from  each  source.  The  implemented  plume  rise 
algorithm  is  an  adaptation  of  Briggs'  formulations  and 
accounts  for  varying  stability,  momentum  and  buoyancy 
rise,  and  stack  downwash.  The  physical  source  param- 
eters required  as  input  to  plume  rise  calculations  in- 
clude emission  type,  source  location,  emission  rate, 
physical  stack  height,  stack  gas  temperature,  stack 
diameter,  stack  gas  exit  velocity,  and  ambient  temper- 
ature and  wind.  A  thorough  description  of  the  plume 
rise  equations  is  available  in  EPA  1980. 

If  appropriate  stack  or  source  parameters  are 
unknown,  a  user-estimated  effective  plume  height  can 
be  entered  directly  into  the  CITPUFF  model. 

Terrain  Height  Adjustment 

Following  the  approach  used  by  other  models  applied 
in  complex  terrain  (Egan  1975,  Schulman  and  Scire 
1980,  Lavery  et  al.  1983),  stability-dependent  plume 
path  coefficients  are  used  to  calculate  the  effects  of  ter- 
rain elevation  on  the  height  of  the  plume  centerline. 

The  terrain-corrected  effective  stack  height,  ht,  above 
the  terrain  is  calculated  from 

h,=  h-(l-T)(minimum{h,ER-Ep}) 

where 

h  =  hs  +  Ah  is  effective  stack  height  above  stack  base 

elevation, 
hs  =  physical  stack  height, 
Ah  -  plume  rise, 
ER  =  receptor  elevation, 
Ep  =  stack  base  elevation, 


22.5°  Sector  Dispersion 

An  option  also  exists  in  CITPUFF  to  choose  the  22.5° 
sector  dispersion  algorithm,  as  applied  in  the  EPA  COM- 
PLEX I  model  (EPA  1981).  The  puff  is  dispersed  horizon- 
tally within  a  22.5°  arc  in  the  direction  of  the  local  wind, 
dependent  on  distance  from  the  source.  Vertical  disper- 
sion is  calculated  identically  to  the  distance-related  az  of 
the  COMPLEX  I  model.  Because  the  COMPLEX  I  ap- 
proach is  applied  for  calculating  az,  the  user  is  cautioned 
that  the  initial  value  of  az  is  assumed  to  be  1  m.  The  user 
cannot  input  initial  az  values  greater  than  1  m  when  the 
22.5°  sector  option  is  applied. 

The  method  used  to  assign  ground-level  concentration 
values  to  grid  receptors  from  the  calculated  sigma 
values  is  the  same  as  that  used  in  all  other  CITPUFF  op- 
tions. In  addition,  the  user  has  the  option  of  specifying  a 
time  at  which  the  time-dependent  dispersion  coefficients 
are  applied. 


Definition  of  Initial  Sigma  Values 

CITPUFF  allows  for  the  specification  of  initial  sigma 
values,  ayo  and  azo,  at  the  pollutant  source  (s  =  0). 
Area/volume  source  distributions  can  be  specified  by 
appropriate  choices  of  the  initial  sigma  values.  With 
some  approximation,  this  allows  the  simulation  of 
various  types  of  complex  sources  including  forest  fires 
and  surface  mines. 


Effective  Emission  Height 

To  complete  the  solution  for  ground-level  concentra- 
tion for  a  discrete  puff  (equation  [10]),  an  effective  pollu- 
tant source  height  (h)  must  be  calculated  from  the 
physical  parameters  associated  with  the  source  emis- 
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t  =  the  plume  path  coefficient,  whose  value  is  given 
in  the  following  tabulation  for  different  stability 

classes: 


Stability  class 

A 
B 
C 
D 
E 
F 


Plume  path  coefficient,  7 

0.5 
0.5 
0.5 
0.5 
0.3 
0.3 


For  neutral  and  unstable  conditions  the  plume  is  lifted 
one-half  of  the  difference  between  the  elevation  of  the 
receptor  and  the  elevation  of  the  stack  base,  with  the  ad- 
ditional restriction  that  the  plume  always  be  at  least 
half  the  height  above  the  ground  that  it  would  be  with  no 
topography.  Under  stable  conditions  the  one-half  condi- 
tion is  replaced  by  a  one-third  condition. 

The  appropriateness  of  plume  path  corrections  to  all 
situations  is  uncertain  and,  therefore,  the  user  has  the 
option  of  inputting  any  desired  set  of  values  if  the 
default  values  are  considered  inappropriate.  Obviously, 
selection  of  7  =  1  will  yield  a  flat,  non-terrain-lifted 
plume. 

Wind  Power  Law  Adjustment 

CITPUFF  uses  a  power  law  function  of  the  form: 


to  adjust  the  winds  input  from  the  WINDS  model  (at 
z  =  zR)  to  the  height  required  by  the  plume  rise  calculation 
for  puff  trajectory  calculations.  The  default  power  law  ex- 
ponents within  CITPUFF  for  each  stability  category  are: 


Stability  class 

A 
B 
C 
D 
E 
F 


Power  law  exponent 

0.15 
0.15 
0.20 
0.25 
0.40 
0.60 


The  user  has  the  option  to  change  the  power  law  expo- 
nent if  desired. 

Effect  of  Mixing  Depth  and  Puff  Fumigation 

Mixing  depth  (H)  is  taken  into  account  by  using  the 
technique  of  multiple  reflections  of  the  puff  when  signifi- 
cant pollution  reaches  the  specified  level  as  first  estab- 
lished by  Bierly  and  Hewson  (1962).  If  the  mixing  depth 
height  is  taken  as  H,  this  technique  leads  to  the  following 
equation  for  the  puff  concentration  at  ground  level: 


x(x,y,0): 
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For  the  case  of  uniform  vertical  distribution  within  H  this 
becomes: 
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In  the  current  version  of  CITPUFF,  the  user  specifies  one 
of  two  possible  algorithms  based  on  equations  [12]  or  [13]. 

Puffs  with  an  effective  emission  height  above  the  user- 
defined  mixing  height  will  not  contribute  to  ground-level 
concentrations.  If,  however,  the  mixing  height  at  some 
later  time  period  becomes  greater  than  the  height  of  these 
puffs,  they  then  begin  to  contribute  to  ground-level  con- 
centrations, although  their  effective  release  height  re- 
mains unaltered.  For  simplicity,  CITPUFF  assumes  that 
puffs  residing  below  the  mixing  depth  can  diffuse 
through  the  maximum  depth  H  encountered  by  the  puffs 
at  any  given  time  of  simulation. 


Puff  Sampling  Function 

The  dosage  integral  in  equation  [11]  for  a  single  source 
is  approximated  in  CITPUFF  by  a  summation  of  the  form 


X  +X  -j    N 

XdX=£ 


X'5x 


[14] 


i=i 


where  the  approximations  x  =  ut  and  dx  =  udt  are 
used.  The  summation  process  is  similar  to  approximate 
integration  by  the  trapezoidal  rule,  except  that  concen- 
tration at  the  end  of  a  time  interval  j,  Xj  is  assumed  to 
prevail  over  the  entire  time  interval  (St.  =  &C./11).  The 
summation  in  equation  [14]  is  obtained  by  periodic 
sampling  of  the  puffs.  Depending  upon  the  number  of  in- 
tervals or  sampling  frequency  used,  the  approximation 
may  converge  to  the  integral  at  any  desired  level  of 
accuracy. 

The  sampling  frequency  algorithm  developed  by  Start 
and  Wendell  (1974)  is  used  to  evaluate  the  distance  in- 
crement form  used  in  equation  [14].  During  periods  of 
strong  winds,  this  algorithm  allows  the  puffs  to  be 
sampled  several  times  during  an  advection  step.  At  low 
speeds,  puffs  will  be  sampled  less  frequently. 

To  determine  the  contribution  of  each  puff  to  ground- 
level  pollutant  concentration,  the  radius  of  each  puff  is 
determined  by  the  relation: 


R„  =  2a 


where 


v    =  the  minimum  concentration  of  interest,  and 

•^■min 

Xp  =  the  concentration  at  the  puff  center. 

The  ground-level  concentration  is  computed  for  each 
time  step  and  accumulated  over  the  entire  time  of  the 
simulation  for  each  grid  point  thaMies  within  the  radius 
of  influence  of  each  puff.  If  xp  ^  xmin  (as  defined  by  the 
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user),  the  puff  is  considered  to  have  become  so  diluted 
that  it  no  longer  can  be  considered  as  contributing  to  the 
ground-level  concentration.  Under  these  conditions  the 
puff  is  eliminated  from  further  calculations  to  conserve 
computer  requirements.  A  puff  that  is  transported 
beyond  the  analysis  area  grid  is  also  eliminated  from  fur- 
ther calculations.  A  maximum  of  500  puffs  can  be  con- 
sidered simultaneously  in  CITPUFF. 


EXAMPLE  APPLICATION  OF  CITPUFF 
Overview 


An  analysis  area  in  northwestern  Colorado  was 
chosen  to  provide  an  example  application  of  CITPUFF. 
Emission  from  two  hypothetical  oil  shale  sources  within 
the  analysis  area  (the  Piceance  Basin)  were  studied. 

Figure  2,  a  map  of  the  analysis  area,  includes  the  loca- 
tion of  the  two  sources  and  nearby  landmarks.  A  sum- 
mary of  the  primary  spatial,  temporal,  and  meteor- 
ological boundary  conditions  applied  in  this  example 
application  of  TAPAS  (WINDS  and  CITPUFF  models) 
includes: 

•  Grid  dimensions:  x  =  81;  y  =  61 

•  Grid  spacing:  2-min  grid  intervals 

•  Approximate  distance  between  grid  points: 
x  =  2.86  km;  y  =  3.70  km 

•  Terrain  and  land  cover  at  each  grid  point  (terrain 
extracted  from  data  tapes  using  the  TAPAS  Ter- 
rain Access  Routines  (NOAA/EDIS/NSGDC  1980); 
land  cover  was  manually  digitized) 

•  Surface  temperature  at  1800  m:  20°  C  (293°  K) 

•  Influencing  wind:  west-southwest  at  4  m/s  (10  m 
above  surface) 

•  Sampling  period:  12-hr  steady-state  wind  condi- 
tion 

•  Assumed  mixing  depth:  4300  m  above  mean  sea 
level 

•  Assumed  stability  for  Pasquill  stability  related 
dispersion:  Class  E  (stable)  (Pasquill-Gifford 
dispersion  coefficients  assumed) 

•  Puff  release  frequency:  10  min 

•  No  deposition  or  decay  of  pollutants  was  con- 
sidered in  the  CITPUFF  model  run 

•  Number  of  minutes  between  puff  sampling:  Deter- 
mined within  model 

•  Ground-level  concentrations  were  output  only  at 
the  end  of  the  12-hr  period  (PTIME(1)  =  12) 

•  Default  values  for  terrain  adjustment  factors  and 
wind  power  law  were  assumed 

•  Puff  dispersion  was  switched  from  distance-  to 
time-related  calculations  at  3-hr 

•  A  12-hr  simulation  with  no  variation  in 
meteorological  conditions  throughout  the  period 
was  performed 

An  elevation  contour  map  of  the  analysis  area  is  pro- 
vided as  figure  3.  A  velocity  vector  plot  of  the  WINDS 
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Figure  2.— Base  map  of  the  analysis  area  considered  in  the  example 
CITPUFF  analysis. 
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Figure  3.— Elevation  contour  map  of  analysis  area  (contour 
interval— 100  m). 

model  results  for  a  west-southwest  influencing  wind 
condition  are  provided  in  figure  4. 

Elevated  sulfur  dioxide  (S02)  emission  data  from  the 
two  oil  shale  sources  considered  in  this  example  are 
listed  below.  To  account  for  the  horizonal  and  vertical 
distribution  of  emissions  from  individual  plant  proc- 
esses, elevated  emissions  from  these  oil  shale  sources 
were  assigned  initial  a  and  ctz  of  125  m  and  75  m,  respec- 
tively. Since  no  engineering  information  was  available 
on  stack  parameters,  CITPUFF  was  not  asked  to  cal- 
culate plume  rise  in  this  example.  All  elevated  emissions 
were  assumed  to  be  emitted  at  a  300-m  effective  plume 
height. 
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Figure  4.— Velocity  vector  plot  of  west-southwest  wind  field  for  the 
example  analysis  area  (based  on  influencing  wind  of  WSW  at  4  m/s). 

Source  V.  Prototype  1  —  101  g  SOz/s  (Cb-MIS  process) 
Source  2:  Prototype  2  —  67  g  S02/s 

(Union-B  DM/SR  process) 

A  puff  trajectory  plot  of  the  analysis  results  is  provided 
as  figure  5.  A  contour  plot  of  ground-level  S02  pollutant 
concentrations  for  the  example  12-hr  period  is  provided 
as  figure  6. 


CONSIDERATIONS  OF  MODEL  VALIDITY 

The  CITPUFF  model  is  designed  to  add  two  features  of 
dispersion  in  complex  topography  that  are  not  incor- 
porated in  other  available  Gaussian-based  models. 
These  are,  first,  the  accommodation  of  a  space-  and 
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Figure  5.— Puff  trajectory  plot  of  the  two  sources  modeled  in  this 
example  CITPUFF  analysis. 
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Figure  6.— S02   ground-level   pollutant  concentrations  (^g/m3) 
resulting  from  this  example  CITPUFF  analysis. 

time-varying  wind  field  and,  second,  the  availability  of 
added  correction  factors  to  adjust  the  plume  path 
beyond  the  capacity  of  the  flow  field.  Both  of  these 
features  are  needed  to  simulate  conditions  in  complex 
terrain.  A  choice  of  three  diffusion  parameterizations  is 
also  provided.  Adjustments  to  diffusion  coefficients  can 
be  easily  accomplished  if  a  consensus  is  established  on 
the  most  appropriate  coefficients.  However,  this  model 
can  not  be  expected  to  resolve  all  controversies  asso- 
ciated with  modeling  dispersion  in  complex  topography. 
This  section  describes  the  capabilities  and  limitations  of 
the  CITPUFF  model.  It  does  not  provide  assurance  that 
the  model  will  always  be  correct  or  that  it  has  been 
totally  validated.  However,  to  the  extent  that  the  model 
incorporates  conventional  procedures  (Gaussian  form, 
Pasquill-Gifford  diffusion  and  its  variants),  one  can 
estimate  its  likely  validity. 

Uncertainty  about  model  performance  should  be  a 
keynote  in  most  model  applications  and  will  remain  so  in 
air  quality  applications  for  the  foreseeable  future  (Fox 
1984).  For  the  particular  case  of  models  that  treat  dis- 
persion for  long  distances  in  complex  terrain,  there  is  no 
adequate  data  base  upon  which  to  evaluate  perfor- 
mance. However,  this  section  will  summarize  the  status 
of  knowledge,  identify  where  CITPUFF  resides  within 
that  status,  and  suggest  some  limitations  on  its  use  and 
the  utility  of  model  output. 


Dispersion  Accuracy 

The  Pasquill-Gifford-Turner  dispersion  formulation, 
upon  which  CITPUFF  dispersion  is  based,  has  been 
repeatedly  defined  as  state  of  the  art  (Pasquill  1974, 
EPA  1978,  Bass  et  al.  1979).  However,  recent  advances 
in  the  understanding  of  boundary  layer  physics 
(Nieuwstadt  and  Van  Dop  1982)  suggest  that  under  con- 
vective  (unstable)  conditions,  the  formulation  is  fun- 


10 


damentally  incorrect,  and  under  stable  conditions  the 
physics  of  the  boundary  layer  may  be  so  complex  that 
simple  models  are  inadequate.  Recently,  Weil  and 
Brower  (1982)  have  incorporated  some  of  these  modern 
concepts  into  a  model.  Model  results  show  significantly 
improved  performance.  A  major  limitation  of  the  Weil 
and  Brower  (1982)  approach,  however,  is  the  need  for 
significantly  increased  meteorological  data  to  drive  the 
model.  The  current  EPA  complex-terrain  model  develop- 
ment program  (Lavery  et  al.  1982,  1983;  Strimaitis  et  al. 
1983)  is  systematically  evaluating  the  effects  of 
topography  on  ground-level  concentration  of  plumes. 
Results  from  this  study  suggest  that  improvement  in  the 
accuracy  of  ground-level  predictions  can  result  from  the 
addition  of  plume  path  correction  factors  and  deter- 
mination of  dispersion  coefficients  based  on  local 
measurements  of  the  turbulence.  Because  local  meas- 
urements are  generally  not  available,  it  is  necessary  to 
rely  on  conventional  dispersion  schemes.  After  consid- 
ering comparisons  of  existing  techniques,  this  study  con- 
cludes that  COMPLEX  I  does  a  better  job  of  simulating 
the  field  data  than  any  of  the  other  routinely  available 
models  (Lavery  et  al.  1982).  COMPLEX  I  is  distinguished 
by  the  fact  that  it  allows  the  plume  to  disperse  horizon- 
tally through  an  arc  of  22.5°  under  all  stabilities.  That 
such  a  simple  approximation  works  is  likely  due  to  the 
presence  of  rather  large-scale  meanders  in  the  flow  in 
complex  terrain  under  all  stabilities. 

The  CITPUFF  model  provides  the  user  with  a  choice  of 
three  available  dispersion  schemes:  Pasquill-Gifford, 
PPSP,  and  22.5°  sector.  Since  all  of  the  schemes  have 
had  extensive  evaluation,  it  is  not  necessary  to  reestab- 
lish their  validity.  Of  concern,  however,  is  how  well  the 
discrete  puff  algorithm  in  CITPUFF  is  able  to  simulate 
the  Gaussian  plume  and,  therefore,  how  well  the  disper- 
sion algorithm  is  implemented.  For  example,  an  impor- 
tant requirement  of  CITPUFF  is  that  it  provide  an  ade- 
quate approximation  of  the  continuous  point  source 
solution  in  situations  where  use  of  the  continuous  point 
source  solution  is  justified. 

Adequate  representation  by  a  puff  model  of  the  con- 
tinuous point  source  solution,  particularly  within  10  km 
of  the  source,  depends  on  the  puff  release  and  puff 
sampling  rates  used.  Guidance  in  selecting  values  of 
these  parameters  has  been  provided  by  an  investigation 
using  the  MESODIF  model  conducted  by  Start  and 
Wendell  (1974).  This  investigation  showed  that  if  the 
puff  travel  distance  between  samples  is  restricted  to 
less  than  2.4  km,  and  the  puff  release  interval  is  no  more 
than  10  min,  the  puff  model  will  provide  an  adequate  ap- 
proximation of  the  continuous  point  source  solution  for 
all  but  isolated  cases  characterized  by  very  strong  hour- 
ly averaged  wind  speeds.  A  release  rate  of  six  puffs  per 
hour  was  found  to  provide  a  reasonable  compromise 
between  the  increased  accuracy  in  representing  the 
continuous  point  source  solution  that  could  be  achieved 
by  more  frequent  puff  releases,  and  the  increased  com- 
putational time  and  costs  involved  in  increasing  the  puff 
release  rate.  For  long-term  model  applications  (on  the 
order  of  a  year)  a  rate  of  three  puffs  per  hour  was 
shown  to  provide  a  reasonable  approximation. 


CITPUFF  (using  the  Pasquill-Gifford  dispersion  option) 
and  the  continuous  point  source  equation  were  com- 
pared to  find  the  normalized  ground-level  concentra- 
tions predicted  by  each  approach  over  flat  terrain 
under  various  stability  conditions.  In  all  cases,  CITPUFF 
was  applied  assuming 

puff  release  rate  =  6  puffs  per  hour 
basic  advection  step  =  10  min 
puff  sampling  time  =  2  min. 

A  wind  speed  of  2.78  m/s  was  tested  under  Pasquill 
stability  classes  B,  C,  D,  E,  and  F.  An  additional  wind 
speed  of  11.11  m/s  was  also  tested  under  the  D  stability 
condition.  Comparisons  were  made  up  to  25  km  from  the 
source. 

Figures  7  through  11  present  the  results  of  the  CIT- 
PUFF versus  continuous  point  source  Gaussian  solutions 
for  Pasquill  stability  conditions  B,  C,  D,  E,  and  F,  respec- 
tively. The  comparisons  indicate  that  CITPUFF  provides 
a  very  close  approximation  of  the  continuous  point 
source  solution  under  the  tested  conditions. 
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Figure  7.— Comparison  of  continuous  point  source  solution  and 
CITPUFF  (Pasquill  dispersion  parameters)  solution  over  flat 
terrain;  stability  B— unstable. 
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Figure  8.— Comparison  of  continuous  point  source  solution  and 
CITPUFF  (Pasquill  dispersion  parameters)  solution  over  flat 
terrain;  stability  C— slightly  unstable. 
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Figure  9.— Comparison  of  continuous  point  source  solution  and 
CITPUFF  (Pasquill  dispersion  parameters)  solution  over  flat 
terrain;  stability  0— neutral 
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Figure  10.— Comparison  of  continuous  point  source  solution  and 
CITPUFF  (Pasquill  dispersion  parameters)  solution  over  flat 
terrain;  stability  E— slightly  stable. 
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Figure  11.— Comparison  of  continuous  point  source  solution  and 
CITPUFF  (Pasquill  dispersion  parameters)  solution  over  flat 
terrain;  stability  F— stable. 


Plume  Path  Correction  Factors 

A  primary  focus  of  current  research  is  on  how  a 
plume  behaves  as  it  encounters  simple  topographic 
elements.  This  research  indicates  that  the  first  require- 
ment for  predicting  reasonable  ground-level  concentra- 
tions is  simulation  of  the  correct  plume  trajectory 
(Lavery  et  al.  1982).  The  puff  trajectories  calculated  by 
the  CITPUFF  model  depend  on  the  wind  field  provided 
by  the  WINDS  model  (or  any  other  method  used  to  gen- 
erate a  grid-based  wind  field).  For  applications  with 
WINDS  the  plume  trajectory  is  limited  by  the  two-dimen- 
sional nature  of  the  WINDS  model.  WINDS  predicts 
flows  parallel  to  the  topography  at  the  ground  surface 
and  hence  contains  no  vertical  velocity  independent  of 
the  orographic  lifting.  Thus,  the  flow  field  is  unable  to 
resolve  any  vertical  adjustment  that  the  puff  may  ac- 
tually undergo  as  it  approaches  terrain.  In  an  extensive 
analysis  of  this  problem  Bass  et  al.  (1981)  concluded  that 
making  stability-dependent  plume  height  adjustments  in 
response  to  the  presence  of  terrain,  is  appropriate. 

Correct  vertical  adjustment  of  the  puff,  however,  may 
require  additional  considerations.  EPA  studies  indicate 
that  the  local  Froude  number  (F  =  U/Nh  where 
N  =  [(g/0)-(d0/dz)],/2)  associated  with  the  height  of  the 
local  terrain  element,  determines  whether  the  plume 
goes  toward,  around,  or  over  the  terrain.  In  fact,  models 
have  recently  been  developed  to  make  precise  numer- 
ical evaluations  of  ground-level  concentration  under 
each  of  these  regimes.  However,  the  needed  meteor- 
ological information  is  not  available  in  the  applications 
of  interest  to  the  targeted  users  of  TAPAS.  Furthermore, 
the  concept  of  precise  plume  behavior  has  not  been 
generalized  to  anything  more  complex  than  a  single 
isolated  hill  considering  wind  speed  (not  direction)  and 
temperature  variations.  Continuing  research  will  likely 
yield  greater  understanding  and  better  data.  Options 
provided  in  CITPUFF  provide  the  model  user  with  the 
ability  to  use  or  not  use  plume  path  correction  factors. 
As  more  precise  plume  adjustment  factors  are 
developed,  they  can  be  incorporated  into  CITPUFF. 

Because  the  CITPUFF  model  operates  with  a  modeled 
wind  field,  it  is  appropriate  to  discuss  the  implications  of 
plume  path  correction  factors  within  this  wind  field.  The 
absence  of  vertical  velocity  in  the  model  leads  to  a 
tendency  for  winds  to  flow  around  topography.  The  path 
correction  features  cause  the  puff  centerline  to  decend 
toward  the  topography.  Thus,  their  use  is  a  conservative 
factor  that  tends  to  cause  higher  values  of  ground-level 
concentration  than  would  otherwise  be  calculated. 

Comparison  of  CITPUFF  to  EPA  COMPLEX  I 
and  COMPLEX  II 

Overview 

Currently  (1984)  EPA  does  not  recommend  any  model 
for  application  in  complex  terrain  (EPA  1978).  Further- 
more, they  suggest  that  appropriate  cautions  should  be 
applied  when  using  Gaussian  plume  models  beyond  50 
km.  However,  in  order  for  the  federal  land  manager  to 
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meet  the  requirement  that  air  quality  related  values  of 
Class  I  areas  be  protected  (Fox  et  al.  1982),  model 
results  must  be  used. 

Although  no  model  is  recommended,  the  EPA  does  sug- 
gest that  COMPLEX  I  and  COMPLEX  II  can  be  used  as 
screening  models.  In  particular,  recent  results  from  the 
small  hill  study  (Strimaitis  et  al.  1983)  indicate  that 
COMPLEX  I  comes  closer  to  predicting  maximum 
ground-level  impacts  than  does  COMPLEX  II.  COM- 
PLEX II  appears  to  be  overly  conservative.  A  compar- 
ison of  CITPUFF  against  COMPLEX  I  and  COMPLEX 
II  for  an  idealized  situation  is  provided  below. 

COMPLEX  I  and  COMPLEX  II 

COMPLEX  I  and  COMPLEX  II  are  multiple  point- 
source  diffusion  models  developed  by  the  EPA  for  appli- 
cation in  rural  and  complex  terrain  (Irwin  and  Turner 
1983).  Both  models  have  been  used  for  regulatory  pur- 
poses. The  models  are  best  applied  within  50  km  of  the 
source.  Beyond  10  km,  mesoscale  influences,  primarily 
spatial  wind  variation,  may  affect  model  results,  par- 
ticularly in  complex  terrain.  For  a  user-defined  time 
period,  the  COMPLEX  models  assume  that  a  single, 
fixed,  wind  direction  and  speed  (defined  at  the  source) 
transports  and  disperses  the  model  plumes.  A  series  of 
options  and  variables  are  available  within  the  models  to 
account  for  plume  terrain  interactions,  stability-related 
wind  profile,  plume  rise,  emission  rate,  receptor  loca- 
tions, and  other  parameters.  The  basic  difference  be- 
tween COMPLEX  I  and  COMPLEX  II  is  the  treatment  of 
horizontal  dispersion  (<jj  COMPLEX  I  assumes  that 
pollutants  are  uniformly  distributed  over  a  22.5°  sector. 
COMPLEX  II  assumes  horizontal  dispersion  based  on 
Pasquill-Gifford  estimates  (Turner  1970).  For  example, 
based  on  these  dispersion  considerations,  COMPLEX  I 
will  provide  lower  downwind  centerline  ground-level 
concentrations  than  COMPLEX  II  under  stable 
conditions. 


Model  Comparisons  on  Theoretical  Flat 
and  Mounded  Terrain 

COMPLEX  I,  COMPLEX  II,  and  CITPUFF  were  all  ap- 
plied to  simulate  ground-level  concentrations  for 
theoretical  flat  and  mounded  terrain  fields.  A  single- 
source  emission  of  100  g/s  with  an  effective  stack  height 
of  300  m  (final  plume  rise)  was  assumed.  An  influencing 
wind  of  4  m/s  from  the  west  (measured  10  m  above  the 
surface)  was  chosen.  A  stable  condition  (Class  E  stabili- 
ty) was  assumed  for  the  comparison. 

Both  flat  and  mounded  simulated  terrain  data  were 
compiled  for  the  test.  The  mounded  terrain  field  con- 
sisted of  two  Gaussian-shaped  hills,  one  centered  at 
30  km  from  the  source  and  the  other  centered  at  69  km. 
A  cross-section  schematic  of  the  mounded  terrain  is  pro- 
vided in  figure  12.  For  CITPUFF  modeling,  the  terrain 
was  defined  within  a  two-dimensional  grid  of  x  =  75, 
y=21,  with  a  distance  between  grid  points  of  1.5  km. 
Receptors  for  COMPLEX  I  and  COMPLEX  II  were  lo- 
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Figure  12.— Cross  section  of  mounded  terrain  data  (or  EPA 
COMPLEX/CITPUFF  comparisons. 

cated  every  3  km  downwind  to  99  km  from  the  source. 

COMPLEX  I  and  COMPLEX  II  were  applied  as  recom- 
mended for  regulatory  purposes.  The  suggested  con- 
figuration of  COMPLEX  I  and  II  (to  conform  to  the  EPA 
complex  terrain  concepts  determined  at  the  Chicago 
Workshop  on  Air  Quality  Models  in  1980,  EPA  1981)  in- 
cludes the  following: 

1.  terrain  adjustment  factors  (0.5,  0.5,  0.5,  0.5,  0.0, 
0.0,  for  A  through  F  stabilities), 

2.  buoyancy-induced  dispersion, 

3.  plume  heights  no  closer  than  zmin  (10  m)  to  recep- 
tors, 

4.  wind  power-law  coefficients  used  in  lieu  of  other 
justifiable  site-specific  data  (0.1,  0.15,  0.2,  0.25, 
0.3,  0.3,  for  A  through  F  stabilities). 

Three  configurations  of  CITPUFF  as  defined  by  varia- 
tions in  the  dispersion  coefficients  (Pasquill-Gifford, 
PPSP,  and  22.5°)  were  compared.  Point-source  emis- 
sions are  modeled  using  initial  ay  and  az  values  of  1.0  m 
and  1.0  m,  respectively. 

Results  of  flat  terrain  comparisons. — The  results  of 
CITPUFF-COMPLEX  centerline  comparisons  for  flat  ter- 
rain are  plotted  in  figure  13.  CITPUFF  (22.5°  sector) 
yields  the  lowest  ground-level  estimates  of  all  modeling 
at  all  receptors.  COMPLEX  I  yields  the  next  lowest 
estimates  of  ground-level  concentrations.  COMPLEX  II 
yields  ground-level  concentration  results  that  are 
generally  higher  by  a  factor  of  five  than  COMPLEX  I 
estimates.  CITPUFF  (Pasquill-Gifford)  estimates  general- 
ly fall  between  COMPLEX  I  and  COMPLEX  II  values. 
CITPUFF  (PPSP)  estimates  are  high  near  the  source,  but 
beyond  66  km  CITPUFF  (PPSP)  estimates  fall  below 
COMPLEX  II  estimates.  Beyond  50  km  from  the  source, 
all  results  presented  fall  within  a  factor  of  five  range. 

Results  of  mounded-terrain  comparisons.— Model  re- 
sults for  mounded-terrain  runs  are  plotted  in  figure  14. 
Basic  differences  exist  in  the  way  COMPLEX  I  and  II 
and  CITPUFF  account  for  terrain  under  stable  condi- 
tions. COMPLEX  I  and  II,  as  applied,  assume  terrain  ad- 
justment coefficients  of  0.0  under  E  stability.  Based  on 
this  assumption,  ground-level  concentrations  for  terrain 
well  above  the  plume  centerline  are  not  calculated.  The 
highest  ground-level  concentrations  are  calculated 
where  the  plume  impacts  the  terrain.  Note  that  no 
values  are  plotted  for  the  COMPLEX  models  at  the  high- 
est terrain  receptor  in  figure  14,  because  the  model  does 
not  calculate  a  value  at  this  point.  To  continue  plume 
transport  beyond  high  terrain,  the  COMPLEX  models 
assume  that  no  terrain  exists,  and  they  continue 
transport  and  dispersion  of  the  plume  effectively 
through  the  mountain  and  beyond  assuming  the  initially 
defined  wind  speed  and  direction. 
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The  CITPUFF  model  depends  on  terrain-related  wind 
fields  (such  as  provided  by  the  WINDS  model)  to  provide 
an  estimate  of  transport  winds.  Variations  in  the  trans- 
port wind  vector  occur  as  terrain  is  encountered.  In  ad- 
dition, a  terrain  adjustment  factor  of  0.3  is  commonly 
used  for  E  stability  in  the  Pasquill-Gifford-Turner  and 
PPSP  versions  of  CITPUFF.  This  factor  indicates  that  the 
plume  is  lifted  approximately  one-third  of  the  difference 
between  the  elevation  of  the  receptor  and  the  elevation 
of  the  stack  base,  with  the  additional  restriction  that  the 
plume  always  be  at  least  half  the  height  above  the 
ground  that  it  would  be  with  no  topography.  Ground- 
level  concentrations  are  therefore  provided  for  all  grid 
points  in  the  plume's  influence.  The  terrain  adjustment 
and  wind  power  law  adjustment  factors  that  are  nor- 
mally used  in  the  22.5°  sector  version  of  CITPUFF  are 
the  same  as  applied  in  COMPLEX  I. 

The  CITPUFF  (22.5°  sector)  model  yields  the  lowest 
ground-level  concentration  in  the  mounded  terrain  test. 
COMPLEX  I  yields  the  next  lowest  values.  COMPLEX  II 
yields  results  a  factor  of  five  or  higher  than  COMPLEX  I. 
CITPUFF  (Pasquill-Gifford)  results  are  similar  to,  yet 
slightly  lower  than,  COMPLEX  II  estimates.  CITPUFF 
(PPSP)  results  yield  high  near-source  values,  but  lower 
values  than  COMPLEX  II  and  CITPUFF  (Pasquill- 
Gifford)  at  distances  beyond  the  second  mound. 
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Figure  13.— Flat  terrain  comparison  of  COMPLEX  I,  COMPLEX  II, 
and  CITPUFF  (plume  height  =  300  m,  emissions  =  100  gls, 
wind— west  at  4  m/s  measured  at  10  m,  stability  =  E,  CITPUFF 
runs  vary  by  the  definition  of  initial  ay  and  az  at  the  source). 
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Figure  14.— Mounded  terrain  comparison  of  COMPLEX  I,  COMPLEX 
II,  and  CITPUFF  (plume  height  =  300  m,  emissions  =  100  g/s, 
wind— west  at  4  m/s  measured  at  10  m,  stability  =  E,  CITPUFF 
runs  vary  by  the  definition  of  initial  <jy  and  az  at  the  source). 


The  location  of  the  highest  ground-level  concentration 
near  the  peaks  occurs  on  the  windward  side  of  the 
peaks.  When  initial  ay  and  az  values  in  CITPUFF  are  set  to 
larger  values  to  simulate  area  sources,  the  location  of 
the  highest  concentrations  could  shift  to  the  lee  side  of 
the  peaks.  These  variations  are  due  to  the  interaction  of 
terrain,  transport  wind,  and  dispersion  parameters. 

Summary 

COMPLEX  I,  COMPLEX  II,  and  CITPUFF  were  com- 
pared over  flat  and  idealized  mounded  terrain  under 
stable  conditions.  Results  were  evaluated  to  a  distance 
of  100  km  from  an  elevated  point  source.  Beyond  50  km 
in  both  terrain  tests,  all  model  results  fell  within  the 
range  of  centerline  ground-level  concentrations  defined 
by  the  low  values  estimated  by  CITPUFF  (22.5°  sector) 
and  high  values  estimated  by  COMPLEX  II.  Results  from 
a  recent  small-scale  field  experiment  (Lavery  et  al.  1982) 
suggest  that  COMPLEX  I  ground-level  values  may  be 
more  realistic  than  COMPLEX  II  values.  CITPUFF, 
depending  on  the  selected  dispersion  option,  estimated 
both  higher  and  lower  values  than  COMPLEX  I  and 
COMPLEX  II. 

Significant  difficulties  in  the  treatment  of  plume  inter- 
actions with  elevated  terrain  also  will  lead  to  differ- 
ences with  model  results.  The  maximum  value  predicted 
by  CITPUFF  for  elevated  terrain  is  also  quite  sensitive  to 
the  assigned  values  of  initial  dispersion  coefficients.  The 
sensitivity  of  near-source  ground-level  concentrations 
due  to  initial  dispersion  coefficient  values  and  due  to 
variations  among  model  concepts  and  applications  are 
well  known  (EPA  1978). 


Validation  Requirements  and  Limitations  on  the 
Applications  of  CITPUFF 

Dispersion  formulations  have  been  extensively  eval- 
uated, and  the  uncertainty  associated  with  their  use  will 
probably  never  be  less  than  a  factor  of  two  at  best  (Fox 
1984,  Fox  et  al.  1983).  In  fact,  this  is  only  under  the  most 
ideal  of  circumstances  very  close  to  the  source.  At 
distances  beyond  a  few  kilometers,  recent  model  evalua- 
tions under  such  ideal  circumstances  (Londergan  et  al. 
1982,  Weber  et  al.  1982,  Bowne  et  al.  1983)  indicated 
that  models  showed  little  skill  in  predicting  an  observa- 
tion at  a  particular  time  and  location.  Performance 
measures  (Fox  1982)  that  reduce  demands  on  the  models 
indicate  that  the  magnitude  and  frequency  of  occur- 
rence of  maximum  values  can  be  predicted  with  some 
skill. 

In  complex  topography,  air  quality  evaluations  have 
tended  to  focus  on  the  near-source  problem,  e.g.,  within 
a  2-5  km  range.  This  is  highlighted,  for  example,  by  the 
fact  that  the  two  data  bases  EPA  currently  recommends 
to  evaluate  complex  topography  models  have  no  moni- 
tors farther  than  5  km  from  the  source.  The  CITPUFF 
model  is  not  designed,  and  should  not  be  used,  for  such  a 
close-in  impact  analysis.  At  these  distances,  the  puff  ap- 
proximation to  the  plume  is  questionable;  the  wind  field 
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is  not  properly  resolved  and  the  puff  centerline  path  cor- 
rection factors  would  need  to  be  evaluated  with  local 
data  on  Froude  numbers. 

All  indications,  however,  suggest  that  at  further  dis- 
tances from  the  source  (say  beyond  10  km)  the  puff  dis- 
persion within  a  wind  field  approach  of  CITPUFF  is 
needed  in  complex  topography  (Egan  1984).4  A  compar- 
ison between  CITPUFF,  COMPLEX  I,  and  COMPLEX 
II  suggests  that  CITPUFF,  applied  using  Pasquill-Gifford 
dispersion  parameters,  yields  higher  ground-level  con- 
centrations than  COMPLEX  I  but  lower  values  than 
COMPLEX  II. 

These  factors  combine  to  suggest  that  CITPUFF,  in  the 
absence  of  specific  validation  studies,  can  be  expected 
to  overpredict  maximum  ground-level  concentrations. 
Because  of  this,  it  may  be  suitable  for  applications  in  en- 
vironmental impact  statements  and  for  projecting  worst 
case  impacts  in  permit  reviews.  If  applied  to  smoke 
management,  it  should  provide  a  factor  of  safety  by 
overpredicting  likely  impacts.  Considerable  experiences 
will  be  needed  to  determine  if  this  factor  of  safety  is  too 
large  or  not.  CITPUFF  should  be  used  with  extreme  cau- 
tion for  anything  other  than  these  worst  case  screening 
(EPA  1978)  applications. 
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APPENDIX  1 


COMPARISON  OF  VARIOUS  DISPERSION 
SCHEMES  USED  IN  CITPUFF 

Comparisons  among  PPSP  and  Pasquill-Gifford  disper- 
sion techniques  as  applied  in  CITPUFF  were  made  for  a 
(stability  classes  A,  C,  and  F)  and  az  (stability  classes  C 
and  F).  These  results  were  also  compared  to  the  Turner 
Workbook  nomograms.  Sigma  values  were  tested  from  0 
to  100  km  distances  from  the  source.  Figures  Al  through 
A3  present  the  <xy  comparisons  for  A,  C,  and  F  stabilities, 
respectively.  Figures  A4  and  A5  present  the  az  com- 
parisons for  C  and  F  stabilities. 
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comparison  (  a)  for  stability  F. 


comparison  ( a  )  for  stability  C. 


3  5  - 


3.0 


2.5 


E 

£  2.0 

N 

(0 

E 

g> 

CO 

1 .5 


1 .0 


0  5 


Stability  =  F 

Pasquill-Gifford  

PPSP   

Turner  Workbook  □ 


jfe^HEEEEi — i-"-i--t--t-J?-T-i--r-i--r-r-r-i — i — i 


0       10       20      30       40      50      60      70       80      90  100 
Distance  from  source  (km) 

Figure   A5.— PPSP,   Pasquill,   and   Turner   vertical  dispersion 
comparison  (  a  )  for  stability  F. 
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APPENDIX  2 


CITPUFF  USER  INSTRUCTIONS 

CITPUFF  has  been  implemented  as  part  of  the 
Topographic  Air  Resources  Analysis  System  (TAPAS). 
TAPAS  provides  a  modeling  environment  that  supports 
the  efficient  application  of  CITPUFF.  Results  of  other 
TAPAS  models  and  data  preparation  routines  such  as 
the  WINDS  model  and  the  Terrain  Acquisition  Routines 
can  be  used  directly  as  input  files  to  CITPUFF.  In  addi- 
tion, the  results  of  CITPUFF  are  generated  in  a  TAPAS- 
structured  file  format  that  allows  them  to  be  directly  ac- 
cessed by  TAPAS  graphic  plotting  routines.  TAPAS  is 
currently  implemented  on  the  Colorado  State  University 
CDC  Cyber  720  computer. 

The  CITPUFF  package  consists  of  a  collection  of  pro- 
grams, data,  and  procedure  files  that  are  designed  to 
configure  CITPUFF  for  various  grid  dimensions,  set 
model  boundary  conditions,  read  data,  and  output 
results  (Fox  et  al.  1984). 

To  run  the  CITPUFF  package,  the  user  must  have  the 
following  files  on  the  computer  account: 

•  CITPUFF— Fortran  code  of  the  CITPUFF  program 

•  PUFPROC— a  multiple  file  of  procedures  to  support 
CITPUFF  configuration,  input,  run,  and  output 
routines 

•  TAPAS  =  structured  files  of  u  and  v  wind  vector 
components  of  simulated  wind  field  data  to  be  used 
by  CITPUFF  in  calculating  puff  trajectories 

•  A  TAPAS-structured  file  of  gridded  terrain  data 
for  the  desired  analysis  area 

•  A  Control  Parameter  File — a  file  generated  by  the 
user  to  define  the  boundary  conditions  for  CIT- 
PUFF modeling 

•  A  Meteorological  Data  File — a  file  defining  the 
date,  time,  mixing  height,  and  ambient  temper- 
ature for  each  hour  of  simulation 

•  USERNBR — a  file  containing  the  user's  computer 
account  number  and  password 

Without  these  files  the  user  will  be  unable  to  run  the 
CITPUFF  package. 


Configuring  CITPUFF  for  the  Appropriate  Grid  Size 

The  computational  arrays  within  CITPUFF  must  be 
configured  for  each  analysis  area.  A  procedure  file  has 
been  developed  to  make  the  configuration  of  CITPUFF 
a  simple  operation.  The  configuration  procedure  (PUF- 
CHNG)  has  the  following  form: 

-PUFCHNG,PUFPROC,MEX,MEY,ICOLSJROWS, 
REL,RUN,STEP,RATE  =  BKGJOBID  -  DD77 

Specifically,  the  user  should  type  the  following: 

1.  A  dash 

2.  PUFCHNG  (This  specifies  that  the  user  wishes  to 
execute  the  PUFCHNG  procedure  file,  which 
directs  the  reconfiguration  of  CITPUFF) 


3.  PUFPROC  (PUFPROC  is  the  group  of  procedure 
files  that  contain  PUFCHNG  and  other  support 
routines  necessary  for  running  CITPUFF) 

4.  MEX  (Number  of  grid  points  in  the  x  direction  of 
the  u  and  v  wind  components  and  elevation  data 
files) 

5.  MEY  (Number  of  grid  points  in  the  y  direction  of 
the  u  and  v  wind  components  and  elevation  data 
files) 

6.  ICOLS  (Number  of  grid  points  in  the  x  direction  of 
the  CITPUFF  computational  grid;  usually  the 
same  value  as  MEX) 

7.  JROWS  (Number  of  grid  points  in  the  y  direction 
of  the  CITPUFF  computational  grid;  usually  the 
same  values  as  MEY) 

8.  REL  (Number  of  hours  of  emission  release  from 
defined  sources) 

9.  RUN  (Number  of  sources) 

10.  STEP  (Number  of  time  steps  in  hours.  This  value 
must  be  at  least  equal  to  the  number  of  hours  of 
the  simulation  plus  one.) 

11.  RATE  =  BKG  (This  space  is  reserved  for  the  CSU 
Computer  Center  job  priority  processing  rate.) 
The  default  rate  is  background  (BKG);  the  user 
accepts  the  default  by  typing  two  commas  („).  If 
other  rates  are  desired,  the  user  must  type: 

PRIME       (Top  priority— full  price) 

BKG  (30%  discount) 

3RD  (Overnight— 40%  discount) 

12.  JOBID  =  DD77  (The  default  JOBID  is  DD77,  in- 
dicating that  the  hard  copy  output  will  be  titled 
with  DD77  followed  by  a  computer-generated 
three-character  code.  The  user  can  accept  the 
default  by  typing  two  commas  („)  or  input 
another  JOBID.) 

When  the  procedure  execution  instruction  is  properly 
typed  on  the  interactive  terminal,  the  user  can  execute 
the  procedure  by  hitting  the  carriage  return.  A  message, 
the  system-defined  JOBID,  will  appear  on  the  screen; 
this  identifies  that  the  procedure  was  successfully  sub- 
mitted. The  result  of  the  PUFCHNG  procedure  will  be  a 
properly  configured  binary  (object  code)  version  of  CIT- 
PUFF, which  will  appear  on  the  user's  account  under 
the  name  PUFFOBJ. 


Control  Parameter  File 

The  Control  Parameter  File  is  created  by  the  user  to 
direct  the  CITPUFF  model.  The  CITPUFF  Control 


19 


Parameter  File  contains  12  card  image  sections  that 
identify  specific  CITPUFF  dimensions,  parameters,  and 
options.  An  example  control  parameter  file  is  provided 
in  figure  A6.  Each  card  type  section  of  the  Control 
Parameter  File  is  listed  in  table  Al  and  highlighted  in 
figure  A6. 


Card 
Type 

1 

2 
3 
4 
5 
6 
7 
8 
8a 


CITPUFF  TEST  DATA 
l.E-10 

1  25  0  7  1  10  2  1  400  1  0  10 
F  F  06  12  18  24  24  24 
2.859000  3.706000  1  81  1  61 

1  81  1  61  1.0  1.0 

2  24  1  "PASQGIF"  3.0 

T 

50.  29.00  27.50  1.0  1.0  5.5  383.0  4.0 
20.  28.75  30.50  1.0  1.0  30.0  325.0  10.0 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 
100.  63. 

T 

.5  .5  .5  .5  0.0  0.0 
T 

.15  .20  .25  .30  .35  .40 

3 

"COPWSW4",  1,  0 
"COPW4",l,  6 
"COPNW3",l,  18 

Figure  A6.— Example  CITPUFF  control  parameter  file. 

Table  A1  —  Card  sections  for  Control  Parameter  File. 


10 

10a 

11 

11a 

12 

12a 


Card  type  Format 


Description 


Free 
(real) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 


TITLE 
CHIMIN 


User-selected  title  of  CITPUFF  model  run 

Minimum  significant  concentration  value  (this  deter- 
mines puff  radius)  (suggested  value  =  1.0E-10) 

Time  in  hours  for  start  of  release 


Total  number  of  hourly  steps  for  simulation  (Set 
NSTEPS  =  NRELSE+  1  in  current  version) 

Decay  rate  constant  for  non-conservative  puffs 
(NC  =  0,  s~ ',  for  conservative  case) 

Month  in  which  release  started 


Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 


Free  KFLSH 
(integer) 


Free 
(integer) 


4  Free 
(logical) 


PRTDOS 


Free 
(logical) 


Free  PTIME(I) 
(integer) 


Free 
(real) 

Free 
(real) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(integer) 

Free 
(real) 


Free  CONVY 
(real) 

Free  NRUNS 
(integer) 

Free  NRELSE 
(integer) 

Free  MNUM 
(integer) 

Character  SIGTYP 


Day  in  which  release  started 


Puff  release  frequency;  number  of  minutes  between 
puffs 

Number  of  minutes  between  each  puff  sampling 
(used  only  when  KVDT  =  0) 

Number  of  advection  steps  into  first  hour  for  release 
(NN 1-1  for  most  cases) 

Maximum  number  of  puffs  carried  before  screening 
out  no  longer  contributing  puffs 

Allows  number  of  minutes  between  puff  sampling 
to  be  variable  as  determined  by  speed  of  puff  and 
specified  minimum  travel  distance  between  sam- 
pling of  puffs.  For  KVDT-1;DT  determined  within 
program.  KVDT  =  0:DT  is  used  as  specified  (sug- 
gested value  =  1) 

Dose  and  concentration  array  controller 

1— initializes  the  dose  and  concentration  array  to 

zero  after  each  print  as  defined  by  PTIME 

0 — allows  the  dose  and  concentration  array  to 

accumulate 

Number  of  minutes  in  basic  advection  step  (sug- 
gested value  =  10) 

Print  Dosage  Flag.  If  set  to  True  (T).  dosage  arrays 
will  be  printed  for  each  source  at  each  of  the 
selected  times  (see  PTIME)  Suggested  value  (F). 
The  total  dosage  at  the  end  of  the  simulation  will 
always  be  printed  regardless  of  the  selection  of 
PRTDOS 

Print  Concentration  Flag.  If  set  to  True  (T),  concen- 
tration array  will  be  printed  for  each  source  at  each 
of  the  selected  times  (see  PTIME).  Suggested  value 
(F).  The  total  concentration  array  at  the  end  of  the 
simulation  will  always  be  printed  regardless  of  the 
selection  of  PRTCON 

I  =  1,6  hours  between  first  six  outputs  of  dose  matrix 
and  concentration  array.  PTIME(6)  is  the  time  be- 
tween all  subsequent  printouts  for  a  given  simula- 
tion 

Spacing  in  kilometers  between  grid  points  in  x 
direction 

Spacing  in  kilometers  between  grid  points  in  y 
direction 

Column  subscript  for  left  boundary  of  concentration 
field 

Column  subscript  for  right  boundary  of  concentra- 
tion field 

Row  subscript  for  bottom  boundary  of  concentra- 
tion field 

Row  subscript  for  top  boundary  of  concentration 
field 

As  for  MSX  on  card  5  but  for  wind  velocity  field 


As  for  MFX  on  card  5  but  for  wind  velocity  field 

As  for  MSY  on  card  5  but  for  wind  velocity  field 

As  for  MFY  on  card  5  but  for  wind  velocity  field 

Conversion  factor  between  concentration  field  and 
velocity  field  grids  (e.g.,  for  CONVX  =  2  grid  spacing 
for  concentration  field  is  twice  that  for  the  velocity 
field) 

As  for  CONVX,  but  can  be  specified  independently 
Number  of  sources  (maximum  50) 
Total  number  of  hours  of  effluent  release 


Model  number:  MNUM  =  1,  Gaussian  distribution  in 
the  vertical;  MNUM  =  2,  uniform  vertical  concentra- 
tion distribution  (suggested  value  =  1) 
Parameter  used  to  select  the  dispersion  method  to 
be  applied:  the  value  for  SIGTYP  must  be  enclosed 
by  double  quotations: 

"PASQGIF"  -  Pasquill-Gifford  based  dispersion  in  y 
andz 

"PPSP"  -  Power  Plant  Siting  Program  method 
"SECTOR"  -  22.5°   sector  dispersion  in  y  and 
Pasquill-Gifford  in  z  (as  applied  in  COMPLEX  I 
model) 
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Free         TMSWTCH  Time  in  hours  at  which  dispersion  calculations  will 

(real)  switch  from  distance-  to  time-based  algorithm.  If 

dispersion  is  to  be  calculated  based  on  distance  for 
the  entire  simulation  or  if  the  SECTOR  dispersion 
option  is  selected,  TIMSWTCH  should  be  set  equal 
to  9999 

Free         CALRISE  Plume  rise  calculation  switch:  T  (true)  Calculate 

(logical)  plume  rise  within  model  based  on  stack  parameters 

given  on  the  source  parameter  cards  (card  type  8a) 

F  (false)  Plume  rise  will  not  be  calculated.  The  final 
plume  rise  is  assumed  to  be  the  term  HS  as  given  in 
the  source  parameter  cards  (card  type  8a) 

Source  Parameter  Cards— One  card  per  each  source  (1  through  NRUNS).  Note:  The 
use  of  the  parameters  listed  below  depend  on  the  selection  of  CALRISE  (card  type  8) 


Free 
(integer) 

Free 
(real) 

Free 
(real) 

Free 
(real) 

Free 
(real) 


HS(NRUNS)  Stack  height  in  meters  above  ground  level 

XSO(NRUNS)  X  coordinate  of  source  in  concentration  grid  units 

YSCHNRUNS)  Y  coordinates  of  source  in  concentration  grid  units 

SIGYI(NRUNS)  Initial  cross-wind  sigma  in  meters 


SIGZCKNRUNS) 


Initial  vertical  sigma  in  meters  Note:  If  the  22.5°  sec- 
tor dispersion  option  is  selected,  the  initial  SIGZO 
must  be  set  as  1,0.  Because  of  the  use  of  the  COM- 
PLEX I  method  of  calculating  o  an  initial  a  greater 
than  1.0  will  not  be  recognizee!  by  the  model.  With 
the  Pasquill-Gifford  or  PPSP  options  any  initial  az  is 
allowed 


The  following  three  values  are  only  used  if  CALRISE  =  T 

EXTVEL(NRUNS)    Stack  gas  exit  velocity  in  meters  per  second 


Free 
(real) 

Free 
(real) 

Free 
(real) 

Free 
(real) 


Free 
(logical) 


EXTTMP(NRUNS)  Stack  gas  temperature  in  "  K 
DIAM(NRUNS)        Stack  diameter  in  meters 


QX(1,NRELSE) 


QX(NRUNS, 
NRELSE) 


Free 
(real) 


Free 
(logical) 


Free 
(real) 


Free 
(integer) 


Emission  rate  in  gas  per  second  (real)for  source 
number  1 

Emission  rate  in  gas  per  second  for  number 
NRELSE 

Terrain  adjustment  factor  (TER)  switch: 

T  (true)  The  user  must  enter  new  terrain  adjustment 
factors  on  the  next  card  for  A-F  stabilities  (card 
type  10a) 

F  (false)  The  default  terrain  adjustment  factors  will 
be  used:  .5,  .5,  .5,  .5,  .3,  .3.  If  NEVVTER  =  F,  card  type 
10a  must  not  be  included  in  the  control  parameter 
deck 

User-selected  terrain  adjustment  factors  for  A-F 
stability  (used  only  if  NEWTER  =  T) 

Wind  power  law  adjustment  factors  (SNP)  switch: 

T  (true)  The  user  wishes  to  select  his  own  wind 
power  law  for  each  stability.  The  user  must  enter  the 
new  power  law  exponent  on  next  card  for  A-F 
stabilities  (card  type  1 1a) 

F  (false)  The  default  power  law  exponents  will  be 
used  for  A-G  stabilities;  .15,  .15,  .20,  .25,  .40,  .60  If 
NEWPOW  =  F  card  type  1 1a  must  not  be  included  in 
the  control  parameter  deck 

User-selected  wind  power  law  exponent  for  A-F 
stability  (used  only  if  NEWPOW  =  T) 

The  number  of  wind  velocity  fields  that  are  to  be 
read  and  used  in  the  simulation 


Wind  Velocity  Field  Cards — One  card  containing  each  of  the  following  parameters 
must  be  included  for  each  wind  velocity  field  to  be  used  during  the  simulation 


Character  FILEAR 
(NUMFIL) 


Free 
(integer) 


LOCARY 
(NUMFIL) 


Free  TIMEAR 
(integer)  (NUMFIL) 


Name  of  the  TAPAS-structured  file  that  contains  the 
velocity  field  u  and  v  arrays  required  for  a  given 
simulation  period.'  File  name  must  be  surrounded 
by  double  quotations 

Block  location  number  of  first  velocity  field  compo- 
nent (u  component)  within  the  file  FILEAR 

Time  (from  start  of  algorithm)  at  which  velocity 
fields  are  to  be  used.  TIMEAR(1)  can  be  zero  or  one 


'A  single  run  ol  the  WINDS  model  produces  a  direct-access  TAPAS-structured  hie  containing 
arrays  for  elevations,  potential  temperature,  divergence,  u  component  of  velocity,  v  component 
ol  velocity,  etc  The  u  and  v  components  are  seguential  so  that  LOCARY  is  set  at  the  location  of 
the  u  component. 


Meteorological  Data  File 

The  Meteorological  Data  File  is  created  by  the  user  to 
represent  all  meteorological  boundary  conditions  except 
wind  speed  and  direction  for  the  simulation.  The  file 
contains  the  date  of  the  simulation,  the  stability 
category,  mixing  depth,  and  ambient  temperature.  The 
file  must  be  at  least  24-hr  long  (contain  24  cards)  even  if 
a  CITPUFF  simulation  is  performed  for  less  than  24  hr.  A 
sample  Meteorological  Data  File  is  provided  in  Figure 
A7.  Table  A2  gives  the  structure  of  the  Meteorological 
Data  File. 
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Figure  A7.— Example  CITPUFF  meterological  data  file. 


Table  A2  — Structure  of  Meteorological  Data  File 


Card  type 

Format 

Parameter 

Description 

1 

Free 
(integer) 

MET(I,1) 

MMDD,  the  first  two  digits  denote  the  month  and 
the  second  two  the  day 

(I  cards  -  one 
card  in  24-hr 
blocks) 

Free 
(integer) 

MET  (1,2) 

Single  digit  from  1  to  6  denoting  stability  class  (1 
unstable,  4  neutral,  6  stable) 

Free 
(integer) 

MET(I,3) 

Height  of  mixing  depth  above  mean  sea  level  in 
meters 

Free 
(integer) 

MET(I,4) 

Average  ambient  temperature  of  analysis  area  in  °K 
(used  only  for  plume  rise  calculations) 

CITPUFF  Model  Run  Procedure  File 

Once  the  CITPUFF  model  has  been  configured,  the 
Control  Parameter  and  Meteorological  Data  files  have 
been  constructed,  and  the  needed  WINDS  model  results 
and  elevation  data  are  present  on  the  user's  account,  the 
CITPUFF  model  can  be  run.  A  procedure  file  has  been 
developed  to  simplify  the  running  of  the  model.  The  pro- 
cedure file  provides  the  interface  between  the  user  and 
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the  CDC  computer  operating  system.  The  run  CITPUFF 
procedure  (RUNPUFF)  has  the  following  form: 

-  RUNPUFF,PUFPROC,CONTROL,METDAT, 
ELEVDAT,CONCEN,PUFFLOC,RATE  =  BKG, 
JOBID  =  DD77.DEST  =  $ID  =  00$. 

Specifically,  the  user  should  type  the  following: 

1.  A  dash 

2.  RUNPUFF  (This  specifies  that  the  user  wishes  to 
execute  the  RUNPUFF  procedure  file  that  directs 
the  execution  of  CITPUFF) 

3.  PUFPROC  (PUFPROC  is  the  group  of  procedure 
files  that  contains  RUNPUFF  and  other  support 
routines  necessary  for  running  CITPUFF) 

The  next  five  parameters  are  specific  to  the  user  ac- 
count and  the  specific  CITPUFF  run.  User  account  file 
names  should  be  substituted  for  the  names  provided  here. 

4.  CONTROL  (The  user  must  enter  the  file  name  of 
the  desired  control  parameter  file) 

5.  METDAT  (The  user  must  enter  the  file  name  of 
the  desired  meteorological  data  file) 

6.  ELEVDAT  (The  user  must  enter  the  file  name  of 
the  TAPAS-structured  file  containing  the  eleva- 
tion data  (in  meters)  of  every  grid  point  in  the 
analysis  area) 

7.  CONCEN  (The  user  must  select  a  file  name  that 
will  contain  the  output  results  of  the  ground-level 
pollutant  concentrations  calculated  by  CITPUFF) 

8.  PUFFLOC  (The  user  must  select  a  file  name  that 
will  contain  the  coordinate  locations  and  other 
puff  identification  information  for  each  puff  from 
each  source) 

9.  RATE  =  BKG  (This  space  is  reserved  for  the  CSU 
Computer  Center  job  priority  processing  rate.)  The 
default  rate  is  background  (BKG);  the  user  accepts 
the  default  by  typing  two  commas  („).  If  other  rates 
are  desired,  the  user  must  type: 

PRIME       (top  priority— full  price) 

BKG  (30%  discount) 

3RD  (Overnight— 40%  discount) 

10.  JOBID  =  DD77  (The  default  JOBID  is  DD77,  in- 
dicating that  the  hard  copy  output  will  be  titled 
with  DD77  followed  by  a  computer-generated 
three-character  code.  The  user  can  accept  the 
default  by  typing  two  commas  (,,)  or  input  another 
JOBID) 

11.  DEST  =  $ID  =  00$  (The  default  destination  code 
for  hard  copy  output  is  the  main  CSU  Computer 
Center  (ID  =  00).  Output  can  be  routed  to  any  CSU 
site  or  specially  designated  location.  To  accept  the 
default  the  user  will  type  two  commas  („).  To 


direct  the  output  to  another  site  the  user  must  in- 
put the  proper  code  as  provided  by  the  CSU  Com- 
puter Center.  For  example,  to  route  all  output  to 
the  RMFRES  the  user  would  enter  $UN  =  AW$. 

When  the  procedure  execution  instruction  is  properly 
typed  on  the  interactive  terminal,  the  user  can  execute  the 
procedure  by  hitting  the  carriage  return. 

CITPUFF  Results 

A  hard  copy  printout  of  all  model  results  will  be 
automatically  directed  to  the  user-selected  destination. 
CITPUFF  will  also  generate  four  new  files  on  the  user's 
account: 

•  PUFF  DAY— This  is  a  dayfile  listing  the  important 
steps  in  the  execution  of  the  CITPUFF  package  and 
includes  total  operating  time  and  cost  related  units. 
The  dayfile  will  indicate  if  any  errors  occurred  in 
the  run. 

•  PUFFOUT— This  is  a  duplicate  computer  file  of  the 
hard  copy  output.  If  a  run  error  is  indicated  in  PUF- 
FDAY,  the  user  should  review  PUFFOUT  to  diag- 
nose where  the  error  occurred  relative  to  the  input 
parameters. 

•  CONCEN— Will  contain  the  concentration  fields  as 
output  in  the  hardcopy  printout  from  CITPUFF. 
CONCEN  (actually  the  user-defined  name  for  CON- 
CEN) will  be  written  to  the  user's  account  as  a 
direct-access,  multiblock,  TAPAS-structured  file. 
Each  concentration  field,  as  requested  by  the  user, 
will  exist  as  a  separate  file  on  the  multiblock  CON- 
CEN file.  A  schematic  of  the  structure  of  CONCEN 
is  provided  as  figure  A8.  CONCEN  can  be  used 
directly  as  input  to  TAPAS  contour  graphics 
packages.  More  often,  CONCEN  is  used  in  conjunc- 
tion with  other  CITPUFF  runs  as  input  to  the  CIT- 
PUFF Results  Data  Handling  Package.  For  a  further 
explanation  of  the  use  of  this  package,  see  the  CIT- 
PUFF Results  Data  Handling  Package  User's 
Guide— A  TAPAS  User's  Guide  available  from  the 
Rocky  Mountain  Station. 

•  PUFFLOC— Will  contain  the  coordinate  location, 
puff  distance,  and  puff  size  of  every  modeled  source. 
PUFFLOC  can  be  used  as  a  direct  input  into  the 
PUFFPLT  graphics  package  which  can  provide  a 
visual  representation  of  the  modeled  plume.  (See  the 
TAPAS— User's  Guide  to  the  Puff  Trajectory 
Graphics  Package  available  from  the  Rocky  Moun- 
tain Station.) 

The  Control  Parameter  File  for  this  example  (COP- 
CONl)  is  listed  in  figure  A9.  The  Meteorological  Data  file 
COPMET  is  similar  to  that  shown  in  Figure  A7  except 
that  the  temperatures  were  fixed  at  293°  for  all  24  hours. 
Wind  field  u  and  v  components  for  this  analysis  were  con- 
tained within  the  TAPAS-structured  file  COPWSW4.  The 
terrain  data  TAPAS-structured  file  was  named 
COPELVC. 


22 


Ground-Level  Concentration  Summary  Files  (in  uglm3} 


Ground-level  concentration  contribution  from  source  1  at  PTIME(l) 


Ground-level  concentration  contribution  from  source  1  at  PTIME(2) 


Source  1 


Ground-level  concentration  contribution  from  source  1  at  PTIME(6) 


Ground-level  concentration  contribution  from  source  2  at  PTIME(l) 


Ground-level  concentration  contribution  from  source  2  at  PTIME(2) 


Source  2 


Ground-level  concentration  contribution  from  source  2  at  PTIME(6) 


Ground-level  concentration  contribution  from  source  n  at  PTIME(l) 


Ground-level  concentration  contribution  from  source  n  at  PTIME(2) 


Source  n 


Ground-level  concentration  contribution  from  source  n  at  PTIME(6) 


All  Sources  for      Total    ground-level    concentration    from    all    sources    for    the  entire 

the  entire  simulation  (NSTEPS) 

simulation 

(NSTEPS) 


Figure  A8.— Schematic  of  CITPUFF  ground-level  concentration  results  file  (CONCERN). 


PROTOTYPE  OIL  SHALE  =  CITPUFF  TEST  DATA 
l.E-10 

1  13  0  7  1  10  2  1  400  1  0  10 
F  F  12  12  12  12  12  12 
2.859000  3.706000  1  81  1  61 

1  81  1  61  1.0  1.0 

2  12  1  "PASQGIF"  3.0 
F 

300.  29.00  27.50  125.  75. 
300.  28.75  30.50  125.  75. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

101.  67. 

F 
F 
1 

"COPWSW4",  1,  0 
Figure  A9.— Control  parameter  file  listing  for  the  example  CITPUFF 

analysis. 


CITPUFF  Configuration  for  the  Example  Analysis 

The  CITPUFF  model  configuration  procedure  instruc- 
tions had  the  form 

-PUFCHNG.PUFPROC.si.ei.si.ei^^s 

The  configured  binary  version  of  the  model  was 
stored  on  the  user  account  as  PFFOBJ. 


CITPUFF  Example  Run 

The  procedure  instruction  to  run  this  example  had  the 
following  form: 

-RUNPUFF,PUFPROC,COPCONl,COPMET, 
COPELVC,CONCENl,PUFLOCl„„ 

The  dayfile  output  from  this  12-hr  simulation  is  listed  in 
figure  A10. 
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11.30.48.  DD77,T5000,PR200. 

11.30.49.  USER.CCC9AA4,. 

11.30.50.  ROUTE,OUTPUT,DEF,DC  =  LP.ID  -  00. 

11.30.51.  ROUTE  COMPLETE. 
11.30.51.RATE.BKG. 
11.30.51.SDAD,  0.001KUNS. 
11.30.51.SDPF,  0.004KUNS. 
11.30.51.SDMS,  0.143KUNS. 
11.30.51.SDCP,  0.014SECS. 
11.30.51.SDSR,  1.036UNTS. 
11.30.51.   RATE  COMPLETE. 

11.30.51.  GET.TAPE6  =  COPCONl. 

11.30.52.  GET.TAPE5  -  COPMET. 

11.30.52.BEGIN,MAKLOC2,PUFPROC,COPELVC,TAPE8. 

11.30.54.IFE,FILE(COPELVC,LO).OR.FILE(COPELVC,PT),LOCAL. 

11.30.54.ENDIF.LOCAL. 

11.30.54.GET.TAPE8  =  COPELVC/NA. 

11.30.56.IFE,FILE(TAPE8,LO),INDIR. 

11.30.56.COMMENT.  *****  GET  SUCCESSFUL.  COPELVC  NOW  LOCAL  AS  TAPE8. 
11.30.56.REVERT.MAKLOC 

11.30.56.  DEFINE,TAPE14  =  CONCEN1.TAPE13  =  PUFLOCl. 

11.30.57.  GET,PUFFOBJ. 

11.30.58.  PUFFOBJ. 
11.35.57.    END  CITPUFF 

11.35.57.    173500  MAXIMUM  EXECUTION  FL. 
11.35.57.    242.599  CP  SECONDS  EXECUTION  TIME. 

11.35.57.  RETURN,  TAPE5,TAPE6,TAPE8,PUFFOBJ,TAPEl3,TAPEl4. 

11.35.58.  REWIND,TAPE10. 
11.35.58.SKIP,EXLAB. 
11.35.58.ENDIF,EXLAB. 

11.35.58.  REPLACE,OUTPUT  =  PUFFOUT/NA. 
11.35.59.STIME    295.667  UNTS. 

11.35.59.  CTIME    245.843  SECS. 
11.35.59.DAYFILE,L  =  PUFFDAY. 

Figure  A10.— Listing  of  the  dayfile  PUFFDAY  for  this  example  CITPUFF  analysis. 


24 


cd 


ed  B 


b 

CD 

I 

Ch 

0) 

a 

x 


CJ  CO 


DO  »  W 


a 

CO 


09 

T) 
id 
-E 

U  <S 


CD 
00 

c 

cd 
_  E 

tH      O  ^ 


0  03 


> 
cd 

-  3 
"2  3 
•S  'a 
?  o 

6 
e 

o 
u 
u 

CO 


2 

"05 


CO 


*s  2  H  1 
-2  ^'3 

u  o  cd 

"3  :>> 

os  ft  H 
03  co  cd 

i.s 

m  _E 

03 


[V]  ~G 


-ED.S 


cd  +j 
-E  1-1 

■4-1  • 

n 

H  X3 


a 

E 

o 
u 

cd  .2 
(-1 

—  oo^_ 
"O  -«  03 

e  a"a 

cd  co  O 

-r"3  fl 


CD 


o  < 

fa  .. 

•  fa 

O  fa 

P  fa 

_•  <D  Ch 

O  0? 


CO 


<tf  B3 

ea- 
rn u 


i-i  in 
.  oo 

CO  05 

fa 


S<  cd 

CO  h— I 
Oh  EG 


co  S 
Ch 

8  9 
ft  "E 

co  cd 

- 1— «  Ch 

TO  ch 

CD 

CD  "V 

g>  CD 

£  "S 

a  o 

Cd  u 

CO  * 

ft  15 

Eh  <» 
1—1  cd 

u| 

u 


cd  = 

S3  3.S  a 

e  -a 
o  g 

o  cd 


E 

...  O 
DO  CD  -Jp 

•S3  2 

rrj  cd  +j 

b  o  e 

-E  ;E  03 

.S  afi 

cd  O 


CD  CD 


T3 

CD 
CO 
B 


cd  co  a  >5> 
SG  -22  3^ 


CD 
CJ 
Ch 

0 
o 

co 

HH  2 

o  E 

>»J2  cd  - 
OB  fa  +3  E 


b  co 
o 

CD 


CO 


CD  DO 

-a  o 


o  a 


c 
3 
o 

a 
a" 

_o 

'■B 

3 

=3 
O 

ft 


DO 

b 
'cd 

-B 
O 


b 
_o 

"co 

tH 

CD 

ft 


CO 

O 


3 


CO  3 

23 


5  u 

Cd  (H 

'  cd 

CD 
CO 
CD 
fa 


£•3 


a 

cd 


C 

CD 
J 

09 

a 

x 

CD  nj 
O 

E03 
DO 
CD  S 

M«2 


<2  3  g 


1°  a 

CJ  Cd 


T3  ffl  u  o 
"cd 

CJ 


CO 


cp  co 


CO    03  CD 

•-•H  tH 

O 
I-h 


•■3  O 


< 
P 

O  c/3 


O 

S 

s 

o 

CJ 

cd  *i 


CO 


M  ft" 

cd  co  cd 
O   O  O 


X 
"ft 

E 
o 

o 


CD 


T3 
C 
cd 

cd" 

CD 


ft  „ 
cd  .2 

M 

D0_ 


u  ^  .5 

•  rH 

S  ft  2 


CD 


5  s 


i  T3 


cd  3 


cd 

co  O 

3  e 

03 


co  X   &  • 

s  J -a  ,2 

ra  ft  o  o 


O  < 
(Jh  .. 

P  Ph 


H  2S 


O  u  -g  06 

<D 
CJ 

C 

o 

CJ 


P  in 
.  co 

CO  Oi 


^  ^H 

CM-- 

rH  U 
CD  ^ 

o 

Ph 

(3 
O 


.  op 

o  S  2 
ft  50  S 
o 

DO  CD  'S3 

■a  -S3* 
3  T3  cd  a 

ft^-2^  g 

C  g--^  cd  o 
3  B  co  o 
Q.  O   03  CD 

cd 


cd  a 
ft -3 

co  cd 

*fH  Ch 
T3  tH 


CD 


CO 


CO 


O 


CD 

Cd  <4-H 

ft  CVD 


Cm 
IJh 

D 

Oh  T? 


CD 

u 

Ch 

B 

o  • 

co 


H   03  03 

°  § 

U  -B  _3 


CD  CO 
co  .03 

'tH 
O 


'B  g> 

iJ  +H 
CJ  © 


CJ 

0 
'cd 

ft  is 


DO 

B 


B 

03 

tH 

°  ft 

^2  SS 

03  DO 

-a  o 

o  ft 

S  2 


B 
B 
O 

E 
a 

o 


o 
ft 

tH 

'3 

oi 
B 

-13 
O 


03 

ft 

co 


CO 

•a 

u 
O 


09 


CO 


co  B 


-a 

CJ 
M 

CO 

CD 


—    CO  (h 


CD 
ft 

—  X 
op  g  M 


C 

CO 


CO 

-a 

CB 

u  — 

CD 

-  E 


03  CD 


J-  SrH 


CD 
DO 

B 

cd 

Oh 

T3 
B 
cd 

-*— ' 

CO 

03 
Ch 
O 

Cm 


CO 


ra  B 

cj  cd 

14-1  CD 

2  « 

•E 


CO 


B 
B 
O 

+->  " 

cd  ft  O  O 

gpsJU 


ft'E 

E  fa 
CO 


E 
E 

CD  O 
co  CJ 

+H  * 

cd  -h 

•5  • 

, — i  03 

n 

«  B 

•2  "5 
co  S 


cd  a 

ft  .3 

co  cd 

Ch 
Ch 

CD 


-a 

03 


o  < 
Ch  .. 

•  Ch 
O  Ch 

dS 

O  U 

P  in 

.  oo 

co  03 
Oh 


tt  B 

cN- 

rH  U 


CO  4-i 


X 
_CD 

"ft 

B  E 
ft  o 

o 


co 


CM  "C 

IS!  ° 
S  Ch 

*  B 

CD  .2 

ft  C/D 


CO 


1/5  cB 
2*  B 

-5  &-cd 

B   M-H  DO 

cj  o  cd 

co    ft  « 

03   CO  cd 

or  ft 

B  s 

o  y  O 

co  ^  U 

03    Oh  „ 

B  cd  .2 

-9  Ch 
-h 

B  "no 

cd  co  O 

03  ft 
52  3  03 

cd  B  h 

-  °  H 

"HI 

.E  2 

ft  M  3 
o 

DO  03  4O 

•S3  g 
b  y  b 

HH  iB  03 

"  ftu 

.3  ft  S 

cd  O 

co  cj 
03  03 

h->  ,B  B 


CD 
CJ 
Ch 
3 

O 


CJ 
CO 

ft 
E 


cd 
B 
CT 

Ch 
'3 

00 

B 


T3 

CD 


03 


CD  CO 
co  03 


?,  -B 


ft  13 

c_    03  . 


Ch  Ch 


U 


CD 

E  - 
&3  2 

03    ft  +3 


O  co 


ft 

CO 
_  Ch 
03  00 


E  2 


o 
ft 

Ch 

•i-H 

cd 

DC 

B 

T3 
T3 
O 


B 
_o 

co 

Ch 
03 
ft 


03 

"O 

M 

o 


0) 


—  -H  CO 
{3    CH  03 

8  51s 


CD 


cd  3 


H 

>  o 


I  5 

O  co 

o  2 

cd  B 


cC  "S 

a-3 

M-H  DO 

o  cd 

CD 
ft  B 

co  cd 
flB  ft 

73  E 

O  o 

3  0 
ft  « 

cd  .S2 


h-j  —  D0_  .tS 
"  H  03  — 


CO 


q  cd  ft 

.O  6 
x  ^  o  J2 

Ch         h  "*  — 

.fa  -S  N  -3 

O  fa   _j  rH*  CJ 

P  ft  -2  2  o 

•  b  22 fa 

O  U  tJ  «  E 

■      .    03    Ch  O 

p  in  y  03  .y 

-  OO  B    ft  cd 

co  C33  O  *  h-< 

2  rH  O  Oh  CO 
O 
Ch 


WJ       K*  S 

Ch    ^  O 

03  B  ft 
ft-B 

co  cd  00 

=3  B  -S 
S-3 


CD 


X^ 
03  U 

ft.E 


ft  13 

1—1  cd 
U-B 


cd  " 

B  4^ 


co  O 

%  E 
ft 

0  P 
-o 

E  co 
cd  g 
O 

__03 

^2 

U  E 

X  CD 

ft  2 
ft  S 

cd  O 
u 

03 

03  -a 

A  B 

S;  « 

03  co 
w  03 

■S  T3 
o 

03  "H* 

E  ^ 

B 


CO 

B 

CT 
Ch 

•rH 

CO 

00 

B 


U  CO 


-2  « 


CJ     CD     ft  H->     Ch  h^ 


T3 

CD 

CO 

B 


E 

CD 
Ch 

°  ft 

CO  2 

03  00 

TJ  O 

O  ft 
O 


SB 
o 

3 

Ch 
O 

03 

H-> 

03 

E 


o 
ft 


00 
E 

TO 

o 


E 

_o 

'co 

Ch 

03 

ft 


CO 

T3 

M 

o 


1935  -  1985 


Rocky 
Mountain  Station 


